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INTRODUCTION 


Many  commonly  occurring  physical  and  chemical  dynamic  systems  have  widely 
separated  time  constants.  These  systems  are  often  represented  by  sets  of 
initial-value  ODE  which  possess  variables  that  rapidly  change  during  time 
intervals  much  smaller  than  the  duration  of  the  phenomenon  of  interest. 

This  presents  the  numerical  integration  difficulties  associated  with  such 
"stiff"  systems.  Thus  even  integration  routines  stable  for  any  step  size 
(so-called  A-stable  methods)  have  accuracy  problems  in  following  the  eigen- 
values large  in  absolute  value  which  damp  out  early  in  the  solution.  These 
errors  can  easily  propagate  to  destroy  the  remainder  of  the  transient. 

The  stability  limitations  involved  with  most  standard  numerical  techniques 
for  an  n-dimensional  system  is  that  they  require  max|hXi| , i = 1,2, ...,n, 
where  the  X^'s  are  the  local  eigenvalues,  to  be  bounded  by  a single  small 
number,  typically  in  the  1 to  10  range.  Thus  if  a single  eigenvalue  is  large 
in  absolute  value,  severe  restrictions  are  placed  on  the  integration  step 
size.  Depending  upon  the  length  of  the  solution  interval  of  interest,  this 
can  demand  a great  deal  of  .computation  time.  Further,  there  are  limits  on 
how  small  h can  be  before  roundoff  errors  accumulate  and  render  the  calcula- 
tion meaningless  [Lapidus  and  Seinfeld,  1971). 


The  practitioner  is  usually  unaware  of  the  nature  of  stiff  systems  and  the 
associated  numerical  integration  difficulties.  Even  arbitrary  application 
of  stiff  methods  is  deemed  significant  enough  in  many  areas  of  application  to 
be  suitable  for  publication.  On  the  other  hand,  typical  solution  character- 
istics and  the  requirement  on  their  elucidation  may  not  be  fully  appreciated 
fay  thti  fiu.r.ir' ~=>i  analyst. 


In  order  to  make  the  nature  of  the  problem  clearer,  SbMsiuei:  ^ ~~“^ific  linear 
time  invariant  evstem 


r _ 

1- 

- 

— 

• 

X 

0 

X 

x (0) 

= 

9 

y 

_K  X2_ 

y 

y (0) 

with  the  analytic  solution 

x(t)  = xtOexpIX^t] 


y (t)  = C^exp [X^t]  t 02exp[X2t] 
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Now  pick  Xj  « X2  < 0 and  Ci  = C2  = 1.  In  this  case  the  contribution  to  the 
solution  of  Xj  is  negligible  after  a very  short  time  period;  yet  its  presence 
will  fix  the  maximum  allowable  step  size  through  the  domain  of  interest  by  the 
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bound  on  max|hAjJ.  This  domain  would  ordinarily  be  determined  by 
i 

mini  A . | 
i 1 

until  a steady-state  has  essentially  been  reached;  for  nonlinear  problems  this 
domain  may  not  be  obvious. 

In  the  present  work  we  have  used  two  completely  different  approaches  for  de- 
veloping feasible  algorithms  for  solving  stiff  differential  equations.  De- 
noting the  stiffness  ratio  (S.R.)  by 

I max ( A . ) I 

Stiffness  Radio  = S.R.  = i 

I min (A . ) I 

I i 1 I 

as  the  ratio  of  the  maximum  to  minimum  eigenvalues,  we  can  define,  in  an  ap- 
proximate way, 

Non-Stiff  Systems 
Moderate-Stiff  Systems 
Strongly-Stif f Systems 

In  the  case  S.R.  < 102  there  are  many  feasible  and  optimal  numerical  algorithms 
in  the  literature  [see,  Byrne  and  Hindmarsh,  1975;  Shampine,  Watts  and  Davenport, 
1976;  and  Enright  and  Hull,  1976];  therefore,  we  shall  not  consider  these  fur- 
ther. When  S.R.  > 1010  the  present  work  has  developed  a singular  perturbation 
technique  which  seems  quite  feasible.  When  102  < S.R.  < 10^°  the  present  work 
has  developed,  in  a preliminary  way,  new  semi-implicit  Runge-Kutta  methods  which 
are  extremely  useful  and  competitive  with  any  other  algorithm  we  have  encoun- 
tered. 


S.R.  < 10 
102  < S.R.  < 1010 


S.R.  > 10 
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I.  A Singular  Perturbation  Approach 


Consider  a two-variable  set  of  first-order  ordinary  differential  equations 
with  a small  parameter  e multiplying  the  derivative  of  one  of  the  variables,  the 
type  of  system  for  which  singular  perturbation  methods  have  been  developed. 


dx 

dt 


f (x,y,e) 


x (0) 


(1) 


g(x,y,e) 


y (0)  = n 


where 


f (0,0, e)  = g (0,0, e)  = 0 

If  (1)  is  linearized  along  its  trajectory,  it  may  be  expressed  as 


1 


=-■  ” 


dx  f 

dt  x 


9x/£  9y/G  y 


Examination  of  the  Jacobian  eigenvalues  indicates  that  their  spread 
increases  the  smaller  e becomes,  one  approaches  zero  while  the  other  grows 
larger  in  absolute  value.  The  occurrence  of  an  eigenvalue  large  in  absolute 
value  defines  the  stiff  problem;  thus  (2)  can  be  regarded  as  the  linearized 
representation  of  a stiff  system  with  widely  separated  eigenvalues.  It  then 
follows  that  such  a stiff  system  and  the  singular  perturbation  form  of  (1) 
are  at  least  locally  equivalent. 

This  equivalence  allows  recent  developments  in  singular  perturbation 
theory  to  be  used  in  obtaining  an  effective  procedure  for  the  numerical  in- 
tegration of  either  equation  type.  It  can  be  shown  that  the  resulting  al- 
gorithm does  not  require  identification  of  the  perturbing  parameter  e,  hence 
is  quite  applicable  to  the  general  stiff  system. 

Consider  initial  value  problems  of  the  form  of  (1)  with  the  perturbing 
parameter  e very  small.  First  assume  the  solution  may  be  approximated  by  the 
simple  first-order  expansion  in  the  (outer)  variables 

x*  = x + ey 

(3) 

y*  = yo  + ZY1 


Substitute  (x*,y*)  into  (1)  for  (x,y)  and  expand  about  (xq (t ) , y0 (t) j • 
ing  terms  with  like  powers  of  e results  in 


Match- 


f(we) 


0 = *(VVe) 


xQ(0)  = ? 


dT  = fx(VY0)xl  + fy(VY0)yl  ' Xl(0)  = ° 


“dt”  = 9x(VY0)xi  + gy(VyO)yi  ' Y0(0)  = n 


!T=°  , yi(o>  = o 

An  inconsistency  can  arise  when  (4b)  is  not  satisfied  by  (C,n).  To  alleviate 
this,  additional  (inner)  variables  are  introduced  which  are  particularly  im- 
portant to  the  very  early  stages  of  the  transient.  Expanding  these  variables 
to  first-order  in  c and  adding  them  to  (3)  gives  the  new  solution  approximation 
(x*,y*) : 


f . ▼ 


' 
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x*  (t ) = x (t ) + ex. (t)  + X (t/e)  + ex, (t/e) 


y*(t)  = yQ(t)  + ey^t)  + Yq  (t/e)  + eY^  (t/e) 

A boundary-layer  type  characteristic  is  imposed  on  the  inner  variables: 


1Ln  xo  • X1  • Yo  ■ Yi  ■ 0 


t/e 


Let  x = t/e  and  make  this  variable  change  in  (1) 

, x (0)  = e 


dX  £ , X 

57  = ef (x,y,e) 


!*Z  = 

dx 


g(x,yfe) 


y(0)  = n 


(6) 


(7) 


(8) 


Now  substitute  (x*,y*)  in  (8)  and  expand  about  (x0(ex)  + X0(x),  y0(ex)  + Y0(x)) 
Matching  terms  with  like  powers  in  e results  in 


^0 

dr 

dY, 


= 0 


(9) 


~dT=  g(xo(eT)  + xo(T)'  yo(£T)  + VT)) 


and 


h 


dx 

~Stm  f(x0(eT)  + xo(T)'  VET)  + Y0(T))  " f(x0(eT)r  y0UT)) 


dY, 


(10) 


dT=  gxXl(T)  + Vi(£T)  + Vl(T)  “ gx{x0'y0)xi(eT) 


- gy(X0'y0)yi(eT) 

Equations  (4)  and  (9)  share  the  initial  conditions 


Note  that  as  a result  of  (7)  and  (9) 


The  conditions  under  which  (6)  may  be  expected  to  be  a valid  solution 
representation  over  the  domain  of  interest  may  be  found  in  Hoppensteadt 
[1971]  . 


Based  upon  these  concepts  a numerical  procedure  has  been  developed  [Aiken 
and  Lapidus,  1974,  and  Miranker,  1973],  which  solves  the  stiff  set  of  equations. 
Since  the  details  are  in  the  literature,  we  present  here  only  a summary  of  the 
most  important  results. 


The  solution 
outer  (xq,  yo)  and 


is  given  in  terms  of  the  zeroth-order  inner  (Xq,Yq)  and 
first-order  inner  and  outer  (Xj , Yj ,x^ ,yj ) terms  (see  (w) ) . 


x = Xx  + exi  + xQ  + exj^ 

Y - yq  + eY1  + yQ  + eYl 


(13) 


where  e is  an  artificial  bookkeeping  indication  of  the  degree  of  stiffness 
defined  by 


ey  = g(x,y,e)  (14) 

where  g = ew.  This  parameter  need  not  actually  exist  or  be  identified.  The 
outer  terms  are  of  more  interest  than  the  inner  terms,  which  are  important 
only  within  a relatively  small  boundary  layer  region  of  the  transient.  For 
systems  stiff  enough  to  require  special  integration  techniques,  the  zeroth- 
order  outer  approximation  often  is  sufficiently  accurate 


x0  = f (x0,yQ)  , x0(0)  = x (0)  (15a) 

o - g<*0»y0>  **  W<X0'Y0>  <15b> 

The  last  equality  in  (15b)  is  made  since  e is  not  zero.  This  is  properly  what 
has  been  referred  to  as  the  pseudo  steady  state  approximation  (pssa) . The 
conditions  for  the  validity  of  (15) , or  for  regular  degeneracy  to  the  low-order 
solution,  briefly,  require  that  the  initial  conditions  x(0)  = C,y(0)  = n be 
within  the  region  of  asymptotic  stability  of 


|^=g(o,y)  (16) 

where  t = t/e,  and  x is  replaced  by  some  constant  vector  a at  each  instant. 

Consideration  of  the  pssa  as  the  zeroth-order  approximation  (5)  reveals 
that  the  region  of  applicability  corresponds  to  the  region  where  the  outer 
variables  are  much  more  dominant  than  the  inner  ones.  The  inner  variables  are 
then  important  only  within  a narrow  initial  boundary  layer  and  thus  can  be 


used  to  define  this  region.  The  zeroth-order  inner  term  for  the  stiff  variable 
is  by  far  the  most  dominant  [Aiken  and  Lapidus,  1974], 


YQ(t)  = Y0(0)exp[9w/3y(x0,y0)t] 


(17) 


where  Yg(0)  = y(0)  - yg  (0) . Experience  by  the  authors  has  indicated  that  (17) 
is  capable  of  providing  an  a priori  estimate  of  the  boundary  layer  for  linear 
and  nonlinear  applications.  In  this  way  the  boundary  layer  is  defined  as  a 
fractional  decay  of  the  zeroth-order  stiff  inner  variable,  the  effective  boun- 
dary layer  length  given  by  t^ 

|*n[Y  (t)/Y  (0)] | 

tT  = — 5 (18) 

ll^xo(0)'yo(0)jlI 

where  | | * | | is  a suitable  matrix  norm.  Since  for  any  matrix  A 
I U|  | > P (A) 

where  p (A)  is  the  spectral  radius  of  A 

p (A)  = max | A^ | 

Vi 


a conservative  estimate  of  the  boundary  layer  is  thus  provided  by  the  use  of  the 
spectral  radius  for  the  matrix  norm,  that  is,  if  the  first  step  can  be  taken 
greater  than  this  boundary  layer  estimate,  the  pssa  is  applicable. 


The  accuracy  of  using  only  the  zeroth-order  approximation  is  indicated  by 
the  magnitude  of  the  first-order  outer  terms  [Aiken  and  Lapidus,  1974] 


ex  (t)  = (ex  (0)  + — ) exp  (at) 

X X cl  cl 

b 


cyi(t)  - t 


yy*’ 

W 


f = 21 

y 3y 


(19) 


etc . 


where 


w f 

„ x y 

a = f - *- 

x w 

y 

w ff 

‘---v 

w 


ex 


• Y (0) 

1(0)  = 2wTU(0V,7T0)]  f(x<0>'Y<°))  - f(*<0),y0<0)) 


(20) 


all  derivatives  are  evaluated  at  (xg,yg),  and  dimensional  notation  has  been  sus- 
pended. Note  that  w may  not  be  singular.  Since  (19)  may  easily  be  evaluated 
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periodically  at  any  time  during  the  solution,  a convenient  upper  bound  on  the 
error  of  using  the  pssa  for  many  common  systems  is  exj (t) /xo (t)  or  eyi (t)/yg (t) . 
If  these  ratios  are  less  than  say  0.001,  the  accuracy  of  the  pssa  is  indicated 
to  be  better  than  0.1%. 

The  preceding  analysis  is  useful  only  for  systems  which  have  w < 0,  an 
initial  monotonically  decreasing  boundary  layer  [see  Aiken  and  Lapiaus,  1975a]. 
Fortunately,  this  seems  to  be  true  for  the  great  majority  of  applications  in 
stiff  chemical  kinetics.  It  also  appears  that  within  this  practical  context, 
the  dependent  variables  divide  into  the  stiff  and  nonstiff  groups,  and  these 
are  often  identifiable  from  a priori  considerations.  If  not,  a few  small  inte- 
gration steps  within  the  boundary  layer  may  reveal  those  variables  with  com- 
paratively rapid  transients. 

A special  characteristic  of  kinetic  systems  is  that  rarely  does  the  model 
represent  the  chemical  phenomena  closely  enough  to  require  better  than  moderate 
accuracy  in  the  numerical  integration.  This  suggests  the  permissibility  of  a 
model  approximation  like  the  pssa.  Thus  for  systems  too  stiff  to  be  integrated 
by  conventional  means,  the  pssa  is  likely  to  yield  quite  adequate  solution  ac- 
curacy. 

The  pssa  may  prove  invaluable  for  the  integration  of  large  systems,  as  ex- 
plicit routines  may  be  used  to  eliminate  the  need  to  invert  a large  Jacobian, 
necessary  in  all  implicit  methods.  When  w is  linear  in  y,  often  the  case  in 
kinetics,  a decomposition  may  be  effected  to  decrease  the  dimensionality. 

Aiken  and  Lapidus  [1975b]  have  also  shown  how  the  initial  conditions  of 
the  specific  system  may  be  chosen  to  eliminate  the  stiff  variables  or  those  with 
large  eigenvalues.  When  numerically  examined  on  a set  of  nonlinear  problems, 
the  strength  of  the  present  algorithm  was  confirmed. 

It  must  be  pointed  out  however  that  the  crucial  point  in  the  use  of  this 
singular  perturbation  approach  lies  in  the  ability  to  decompose  an  initial  set 
of  ODE  into  the  stiff  and  nonstiff  form  of  (1) . When  the  original  system  has 
a large  dimension  (n  > 10) , such  a discrimination  may  not  be  obvious.  Further, 
the  eigenvalues  of  the  original  system  must  cluster  in  groups  rather  than  be 
spread  out  over  roughly  equal  intervals.  When  this  happens  and  the  S.R.  > 106  - 
1010,  the  algorithm  is  an  extremely  efficient  procedure  for  solving  stiff  ODE 
(see  comments  in  later  discussions) . 


II.  Semi-Implicit  Runge-Kutta  Methods 


At  the  same  time  it  must  be  recognized  that  the  above  singular  perturbation 
approach  has  certain  system  restrictions;  thus  there  is  a question  as  to  whether 
it  can  serve  as  the  format  for  a general  purpose  algorithm  for  solving  stiff 
ODE.  As  a consequence,  our  work  has  proceeded  along  what  might  be  called  more 
conventional  directions  but  with  a special  emphasis.  Here  we  present  some  pre- 
liminary results  on  the  development  and  use  of  imbedded  semi-implicit  Runge- 
Kutta  methods  with  special  error  monitoring  characteristics.  For  problems  with 
102  < S.R.  < 1010  this  approach  seems  to  be  the  most  efficient  that  we  have  en- 
countered. 





~ 
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As  pointed  out  by  Lapidus  and  Seinfeld  [1971] , Runge-Kutta  integration 
techniques  may  be  classified  as  explicit  or  semi-implicit  or  implicit;  the 
explicit  and  implicit  forms  may  be  discarded  as  viable  techniques  for  solving 
stiff  equations  either  because  of  extreme  stability  (and  thus  step  size) 
bounds  or  the  high  level  of  iteration  required.  By  contrast,  semi-implicit 
methods  are  A-stable  and  require  no  iteration.  Thus  this  class  of  methods 
remains  as  possible  candidates  for  a general  purpose  algorithm. 

Perhaps  the  best  semi-implicit  Runge-Kutta  algorithm  developed  to  date  is 
due  to  Michelsen  [1976]  which  we  show  below 


?n+l  " ?n  + Rlh  + R2k2  + V?3 


(21) 


with 


k = h [I  - ha  J(y  )]  f (y  ) 
- I.  _n  _ ~n 


k = h [I  - ha  J(y  )]  f (y  + b k ) 
-2  _ 1,  _n  _ _n  2~1 


J (y  ) = Jacobian  matrix 
~ in 

at  y 
~n 


(22) 


*3  - htl  - ♦ b32k2! 


Given  the  solution  yn  at  xn  the  solution  is  advanced  over  the  increment  h 
to  xn+^  to  yield  Yn+i  v^a  • Equation  (21)  can  be  used  once  kj  , £2,  and 
k3  are  calculated  serially  assuming  all  the  parameters  Rj , R2,  R3,  aj,  b2, 
b3i  and  b32  are  determined.  This  latter  feature  can  be  handled  by  matching 
to  Taylor  Series  expansions  and  using  exponential  fitting.  Thus  Michelsen 
determined 


ax  = 0.4358... 


„ 11 

R1  = 27  " b31 


„ 16 
R2  27  " b32 


r3  = ! 


31 


’ b2  ■ ! ’ b32  ■ 9^-  l6al2 

gi-  (8.32  - 2a3  . 1) 


^1  + 1} 


(23) 


These  parameters  make  (21)  and  (22)  A-stable  and  even  further,  strongly  A-stable. 

However,  it  is  necessary  to  add  a step-size  adjustment  feature  so  that 
when  y is  changing  rapidly  h can  be  decreased  and  vice  versa.  Only  with  this 
adjustment  can  the  algorithm  become  truly  efficient.  This  is  usually  done  by 
the  one-step/two-step  extrapolation  in  the  form 


one-step 


n+1 


— 


y 


such  that 


- v(.h'2'  * 


:n+l  ~n+l 


2?  - 1 


(h/2)  (h) 

Yn+i  - yn+i 


where  p is  the  order  of  the  basic  method.  In  the  present  case  p = 3 and 

(h/2)  1 m 

v = v + — T 

in+1  in+1  7 _n+l 


where 


m (h/2)  (h) 

T = y ,,  - y = truncation  error 

~n+l  in+1  in+1 


Using  (26)  as  an  example,  the  step-size  can  be  adjusted  such  that  ||Tn+1|| 
some  error  bound.  However,  the  amount  of  computation  required  to  go  from 
xn  to  xn+2  has  been  increased  by  200%  over  the  non-error  monitoring  case. 

In  the  present  work  we  have  developed  a completely  different  approach  to 
the  error-monitoring  procedure.  Thus  we  define  a new  algorithm 

?„«.  ■ ?„ + Eih + 1271 

where  kj  and  k2  are  identical  to  those  in  (22) . However,  we  relax  the  order  of 
the  method  by  1 (p  = 2)  by  specifying  b2  = 1 - 2aj  and  then  applying  the  Taylor 
Series  expansions  and  exponential  fitting.  The  end  result  is  a second-order 
method,  (27) , imbedded  in  a third-order  method  (21)  , for  which  (27)  can  be  cal- 
culated at  essentially  no  computer  cost  once  (21) -(22)  have  been  evaluated  over 
the  step  h.  This  second-order  method  is  also  A-stable.  Thus  we  calculate  (21)- 
(22)  to  generate  yn+l*  use  (27)  at  almost  no  cost  to  generate  another  yn+1  and 
compare  the  two.  The  number  of  digit  agreement  can  be  used  to  estimate  the 
truncation  error  and  thus  provide  a complete  error  monitoring  procedure. 

To  illustrate  the  results  obtained,  we  select  the  fluidized  bed  system  de- 
tailed by  Luss  and  Amundson  [1968] 


= i-3(y3  - y1)  + 1.04  x io4ky2 


yx(0)  = 759.167 


= 1.88  x 103(y4  - y2(l  + k)) 


y2(0)  = 0 


s — = 1752  - 269y  + 267yi 

at  3 1 


y3(0)  = 600 


df  " °-‘l  + 320y2  " 321y4 


y4(0)  = o.i 


where  k = 0.0006  exp[20.7  - 15000/yj].  This  system  has  a S.R.  = 106,  and  we 
wish  to  integrate  from  t = 0 to  tf  = 500.  The  computing  time  required  by  the 


BMH 


- 10  - 


one-step/two-step  procedure  is  1.76  CPU  units,  while  that  for  the  present 
imbedded  algorithm  is  0.42  CPU  units.  Obviously  the  present  method  is 
considerably  more  efficient  than  the  one-step/two-step  approach.  We  have  also 
tested  the  current  imbedded  algorithm  vs.  essentially  every  other  algorithm 
in  the  literature;  this  for  small/large  dimensional  systems  (n  = 2 to  50) , 
linear/nonlinear  systems,  stiffness  ratios  of  102  S.R.  <_  1016  and  in 
single/double  precision  arithmetic.  The  preliminary  results  indicate  that 
the  current  algorithm  is  more  efficient  than  any  competitive  procedure. 

As  problems  are  considered  where  the  S.R.  approaches  1010-1012,  the  cur- 
rent semi-implicit  method  may  have  diff iculties  unless  sufficient  computer 
precision  is  allowed.  However,  it  is  in  just  this  region  that  the  singular 
perturbation  approach  of  Part  I of  the  report  becomes  quite  efficient.  Thus 
one  could  suggest  that  the  two  algorithms  developed  in  the  present  work  will 
handle  any  set  of  stiff  ODE. 
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SECTION  I 
INTRODUCTION 

The  Air  Force  Weapons  Laboratory  (AFWl)  Is  Investigating  the  I X F Implo- 
sion of  thin  cylindrical  foils  In  short  z-plnch  devices.  These  devices  are  driven 
by  low  Inductance  capacitor  banks  which  are  electrically  connected  to  the  z-plnch 
by  flat  plate  transmission  lines. 

A computational,  code  (PLATE)  was  written  to  evaluate  two  areas  of  concern 
In  the  design  of  these  capacitor  banks.  First,  uniform  current  flow  Into  the 
azimuthal ly  symaetrlc  load  Is  desired  because  asymmetries  may  cause  Irregular  foil 
Implosions.  Second,  a measure  of  the  capacitor  bank’s  effective  Inductance  Is  use 
ful  since  this  Inductance  must  be  kept  to  a minimum  so  that  the  system  can  dis- 
charge quickly.  PLATE  calculates  azimuthal  current  symmetry  In  the  transmission 
lines,  and  It  estimates  the  effective  transmission  line  Inductance.  To  do  this, 
current  Is  constrained  to  flow  from  simulated  capacitors  through  a square  mesh  of 
Inductors,  capacitors,  and  resistors  that  simulate  the  electrical  characteristics 
of  an  actual  transmission  line.  Important  results  are  displayed  In  calcomp  plots. 
The  major  approximation  In  these  calculations  Is  that  the  transmission  plate  sep- 
aration IS  considered  to  be  small  compared  to  the  mesh  size;  otherwise,  mutual  In- 
ductance effects,  which  are  Ignored,  can  become  Important. 

The  two  objectives  of  this  report  are  to  document  the  code  PLATE  and  to 
provide  a users  manual.  To  accomplish  these  objectives,  both  theory  and  applica- 
tion are  discussed.  Sample  Input  Is  provided  In  all  'Instances  where  an  example  Is 
being  discussed.  A listing  and  sample  output  are  provided  In  the  Appendixes. 
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SECTION  II 

^ EXPERIMENTAL  SETUP 

a. 

The  transmission  lines  used  at  the  AFWL  consist  of  two  parallel  aluminum 
plates  which  are  separated  by  mylar.  Generally,  the  top  transmission  plate  Is 
used  to  carry  current  to  the  load,  and  the  bottom  plate  .Is  at  "ground"  and  acts 
as  a return  current  carrier. 

A typical  experimental  apparatus  consists  of  a rectangular  transmission 
line  with  two  capacitor  bank  modules  attached  to  opposing  sides.  A circular  hole 
of  approximately  10  cm  radius  (R1  In  figure  1)  Is  located  In  the  center  of  the 


Figure  1.  Typical  Square  Transmission  Plate  Problem. 


plate.  The  load  Is  bolted  to  the  periphery  of  this  hole.  Figure  1 Is  a schematic 
of  this  apparatus. 

Outside  a Radius  R2,  the  two  plates  are  separated  by  0.15  cm  of  mylar  (fig- 
ure 2).  The  circular  electrodes  have  a radius  R1  and  separation  of  1 cm.  Between 
the  R1  and  R2  radii,  the  plate  separation  varies  with  radius  depending  on  the 
actual  chamber  design.  In  PLATE  calculations  this  region  Is  assumed  to  have  a 


/ 
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Figure  2.  Cross  Section  of  Transmission  Plate. 


constant  plate  separation.  This  value  of  separation  Is  determined  by  Inductively 
matching  the  actual  chamber. 

"Typical"  values  for  the  variables  annotated  In  figures  1 and  2 are: 

AZi  * normal  plate  separation  ■ 0.15  cm 
a z2  ■ Increased  plate  separation  ■ 1.5  cm 
ex  ■ e2  ■ 2.8eQ  (for  mylar) 

Rj  ■ radius  of  the  load  ■ 0.10  m 

R2  ■ outside  radius  of  Increased  plate  separation  ■ 0.14  m 

NOTE:  For  computational  purposes,  Rj  ■ 0.25  m and  R2  ■ 0.40  m.  The  Increased 
radii  allow  the  respective  arcs  to  approximate  circles  when  they  are  superimposed 
on  a square  mesh. 
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SECTION  III 

THEORETICAL  ARRANGEMENT 


ft 

This  section  explains  how  the  experimental  problem  Is  set  up  as  a calcu- 
latlonal  problem  and  how  the  solution  Is  calculated  as  a function  of  time. 

1.  Problem  Set-up 

A symmetric  one-fourth  of  the  experimental  apparatus  discussed  In  section 
II  Is  computed  by  PLATE.  Figure  3 shows  this  symmetrical  quarter  of  the  system. 


The  quarter  plate  Is  divided  Into  a mesh  of  square  elements.  A schematic 
of  a single  element  Is  shown  In  figure  4. 


Figure  4.  One  Square  Element. 
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The  capacitance  of  a parallel  plate  capacitor  Is  given  In  Lorraln  and 
Corson  (reference  1)  as: 


C ■ e e §- 
r o s 


0) 


where  S Is  the  plate  area,  s Is  the  plate  separation,  ef  Is  the  relative  per- 
mittivity of  mylar,  and  eQ  Is  the  permittivity  of  free  space.  Since  e - ereQ 
and  the  plate  area  Is  aXaY, 


r ■ AXaY 
c-  Z~ZT 


(2) 


This  Is  the  form  of  capacitance  used  In  PLATE. 

The  Inductance  In  an  Infinite  parallel  plate  system  (no  fringing)  can  be 
computed  by  comparing  two  mathematical  expressions  for  the  enclosed  magnetic  field 
energy: 


1/2  LI2  ■ magnetic  field  energy 


/*B2 

J ^ 


dV 


(3) 


With  the  assumption  that  displacement  currents  are  zero,  the  Integral  form  of 
Amperes  law  Is: 


f IT*<rr  • u0  f J‘cTa 


(4) 


Figure  4 shows  the  parallel  plate  situation.  The  first  Integration  Is  accomplished 
along  the  dotted  line.  The  second  Is  done  over  the  area  enclosed  by  the  dotted 
line.  The  Integrations  yield: 


(5) 


1.  Lorraln,  P.  and  Corson,  D.  R.,  Electromagnetic  Fields  and  Waves.  San  Francisco: 
W.  H.  Freeman  and  Company,  1970. 
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Substituting  for  B In  equation  (3)  gives 


u0I2 

• -2—  aXaYaZ 
2aY2 


so  that 


L 


UqAXaZ 

— Er- 


ie) 


(7) 


This  Is  the  form  of  self  Inductance  used  In  PLATE. 

In  PLATE,  square  cells  are  used.  Thus,  the  Inductance  Is  the  same  In  the  X 
and  the  Y directions,  and  only  one  Inductance  array  Is  needed.  Because  the  cells 
are  constrained  to  be  square,  the  cell  capacitance  and  Inductance  are: 


C2  - c(aX)2/aZ 

(8) 

AL2  - UqaZ 

(9) 

where  C2  and  AL2  are  the  PLATE  variables  for  capacitance  and  Inductance. 

Furthermore,  since  mutual  Inductances  are  smaller  than  self  Inductances, 
they  are  Ignored  (see  Appendix  C),  and  the  self  Inductance  becomes  the  total  1n« 
ductance  In  each  cell. 

The  resistance  for  one  cell  Is  calculated  next.  Resistance  depends  upon 
skin  depth,  resistivity,  and  material  geometry. 
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Skin  depth  Is  given  by  Slater  and  Frank  (reference  2)  as 


00) 

where  a ■ conductivity  and  u - frequency.  The  resistance  of  a cell  Is  given  by 

R’pir  (in 

where  p ■ 1/c  Is  the  resistivity,  L Is  the  length  of  the  resistor,  and  A Is  the 
cross-sectional  area  of  the  resistor.  For  a cell  of  width  W,  length  L,  and  skin 
depth  6 (the  skin  depth  Is  much  less  than  the  plate  thickness),  the  resistance  Is 

R “ p W (12) 

In  PLATE,  the  cells  are  constrained  to  be  square;  hence 


R 


L 

6 


Substituting  for  6 and  using  p ■ 1/a 


03) 


R ■ T (14) 

Thus,  the  resistance,  R,  Is  Independent  of  all  geometry  factors  when  square 
cells  are  used. 

The  Handbook  of  Chemistry  and  Physics  (reference  3)  gives  the  resistivity  of 
coomerclal  aluminum  as  2.828  x 10"8  ohm-m  at  20°C.  The  frequency  of  the  current 
wave  Is  on  the  order  of  0.25  x 106  Hz.  Thus,  u - 2irF  ■ 1.57  x 106  Hz  and  the  re- 
sistance of  one  cell  Is  1.67  E-04  ohms.  This  value  Is  used  In  PLATE. 

2.  Slater,  J.  C.  and  Frank,  N.  H.,  Electromagnetism.  New  York;  McGraw-Hill  Book 
Company,  Inc.,  1947. 

3.  Hodgman,  C.  0.,  Weast,  R.  C.,  and  Selby,  S.  M.,  Handbook  of  Chemistry  and 
Physics,  Cleveland:  Chemical  Rubber  Publishing  Company,  1960. 
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2.  The  Mesh 

A circuit  schematic  for  a single  cell  Is  shown  In  figure  5.  Current  Is 
allowed  to  flow  In  both  the  X and  Y directions. 


Figure  5.  Circuit  for  One  Mesh  Element. 

Many  of  these  circuits  elements  can  be* fitted  together  to  form  a represen- 
tation of  an  entire  transmission  plate.  Figure  6 shows  a portion  of  such  an  array. 
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In  figure  6.  each  node  Is  connected  to  a capacitor  which  Is  connected  to 
ground  and  to  four  Inductor  and  resistor  pairs. 

The  theoretical  transmission  plate  has  a load  connected  at  a 25  cm  radius 
circle  and  an  annular  ring  at  40  cm  radius  In  which  the  plate  separation  has  been 
Increased  to  1.5  cm  from  0.15  cm.  Capacitance  Is  Inversely  proportional  and  In- 
ductance Is  directly  proportional  to  plate  separation.  aZ.  Consequently,  the 
elements  of  capacitance  and  Inductance  that  are  In  the  annular  ring  where  aZ2  * 
1.5  cm  have  a decreased  capacitance  and  Increased  Inductance  by  the  multiplica- 
tive factor  Lli/blx. 

Other  changes  In  the  plate  separation  may  be  simulated  by  changing  the  ca- 
pacitance and  Inductance  In  suitable  regions  of  the  transmission  plate  mesh. 
Modeling  of  this  kind  can  be  used  to  make  the  simulated  current  flow  more  symmet- 
ric as  will  be  seen  In  section  IV. 

Inside  the  radius  » 0.25  cm,  the  load  Is  simulated  as  a short  circuit. 
This  Is  done  by  making  the  capacitance  arbitrarily  large  Inside  the  radius  Rj. 

The  capacitor  bank  portion  of  figure  1 Is  simulated  by  an  array  of  capaci- 
tors and  associated  Inductors  connected  to  the  side  of  the  transmission  plate 
mesh.  A schematic  of  the  capacitor  bank  connection  Is  shown  In  figure  7. 


Figure  7.  Capacitor  Bank  Circuitry. 
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The  number  of  simulated  capacitors,  NBK,  In  one-half  of  one  module  Is  equal 
to  the  number  of  mesh  points  Inside  the  length  of  that  half  module.  If  CT  and  ALT 
are  the  capacitance  and  Inductance  of  one-half  of  a module,  then  the  capacitance 
and  Inductance  of  each  simulated  capacitor  Is  determined  by 

Cl  ■ CT/NBK  (15) 

and 

AL1  - ALT  * NBK  (16) 

where,  again  NBK  Is  the  number  of  simulated  capacitors. 

The  problem  Is  set  up  by  Initializing  all  of  the  capacitances  and  Induc- 
tances as  discussed  above.  The  voltage  everywhere  Is  set  to  zero  except  In  the 
capacitor  banks  where  It  Is  set  to  VO,  an  Input  variable.  The  currents  are  every- 
where set  to  zero. 

The  times tep,  OELT,  Is  determined  by  the  formula 

DELT  - • V «.2*C2  (17) 

where  C2  and  AL2  are  the  capacitance  and  Inductance  In  the  main  portion  of  the 
transmission  line.  The  problem  Is  now  Initialized. 

3.  Difference  Equations 

The  two  differential  equations  shown  below  are  used  to  derive  the  difference 
equations  used  In  PLATE.  After  this  derivation,  an  overview  of  problem  solution  Is 
presented. 

al/at  - — [ — (18) 


av/at  ■ I/C 


(19) 


where  V ■ voltage  ■ f(x,y,t) 
I • current  ■ f(x,y,t) 
L • Inductance  ■ f(x»y) 

R • resistance  ■ f(x,y) 

C ■ capacitance  ■ f(x,y) 
t • tine 


Letting  !n*l  and  In  be  the  new  and  old  currents,  respectively,  the  difference 
form  of  equal  ton  (18)  Is: 


The  new  current  Is  solved  for  explicitly 


This  equation  Is  used  in  both  the  vertical  and  horizontal  directions  of  the  mesh. 
The  spatial  relationships  of  the  variables  are  shown  in  figure  8.  In  this  figure 
J„  and  Ju  are  the  vertical  and  horizontal  currents. 


Figure  8.  Spatial  Relationships  of  Currents  and  Voltages 


Equation  (19)  Is  used  for  the  derivation  of  the  second  difference  equation. 
Letting  Vn+^  and  Vn  be  the  new  and  old  voltages,  respectively,  the  difference  form 
of  this  equation  Is 
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V„  + At 
n 


(22) 


The  spatial  relationships  of  the  four  currents  to  a node  and  the  voltage  at  that 
node  Is  shown  In  figure  9. 


Figure  9.  Spatial  Relationships  of  Currents  and  Voltage. 


The  variables  used  In  the  FORTRAN  coding  are  somewhat  different  from  those  used 
here.  Appendix  B should  be  consulted  for  a definition  of  all  Important  variables. 
Appendix  E contains  a listing  of  PLATE. 

Initially,  at  t ■ 0,  the  charge  Is  present  In  the  capacitor  bank;  as  time 
progresses,  It  Is  moved  throughout  the  mesh  and  eventually  absorbed  In  the  short 
circuit  load.  The  problem  solution  proceeds  In  four  parts  during  each  step: 

1.  Equation  (21)  Is  used  to  determine  the  new  currents  through  the  Induc- 
tors connecting  the  capacitor  bank  to  the  transmission  line. 

2.  Equation  (22)  Is  used  to  determine  the  voltage  left  on  the  capacitor 
bank  after  part  1. 

3.  Equation  (21)  Is  used  to  determine  the  new  currents  between  all  nodes 
In  the  mesh. 

4.  Equation  (22)  Is  used  to  determine  the  new  voltages  on  all  capacitors 
after  part  3. 

The  above  solution  scheme  appears  to  be  unstable,  but  It  Is  not.  The  two 
differential  equations  are  coupled  and  are  solved  alternately  In  time.  They,  thus, 
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provide  feedback  to  each  other,  and  the  solution  Is  a variation  of  the  leap  frog 
scheme. 

At  various  Intervals  the  solution  Is  Interrupted  to  allow  printing  and/or 
microfilm  plotting  of  the  data. 

Appendix  D contains  the  results  of  two  sample  problems  which  validate  equa 
tlons  (21)  and  (22)  and  the  solution  scheme  In  which  they  are  used. 
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SECTION  IV 

TRANSMISSION  PLATE  DESIGN 

PLATE  determines  the  current  asymmetry  at  a cylindrical  load  which  Is  at* 
tached  In  the  center  of  a parallel  plate  transmission  line.  PLATE  can  be  used  to 
compare  different  transmission  plate  designs.  The  obvious  comparison  Is  current 
symmetry.  Another  Important  comparison  that  can  be  made  Is  system  Inductance. 

Greater  transmission  plate  Inductance  smooths  the  current  and  produces 
better  current  symmetry.  If  a very  low  Inductance,  high  energy  system  Is  desired, 
transmission  plate  design  may  Involve  a trade-off  between  raising  the  Inductance 
and  Improving  the  current  symmetry.  The  example  explained  In  this  section  In- 
volves exactly  this  situation. 

The  objective  Is  to  compute  the  current  symmetry  of  a 1.2  meter  square 
transmission  plate  with  a load  radius  of  0.14  meter  as  shown  In  figure  10.  There 


Figure  10.  Plate  Schematic  I. 

are  two  capacitor  modules  of  half-width  0.16  meter,  half  capacitance  5.55  tiF,  and 
half  Inductance  of  24.0  nH  (12.0  nH  each  module  because  the  halves  are  In  parallel) 
centered  on  opposing  sides  of  the  plate. 

A symmetric  one-fourth  of  this  system  was  computed  using  the  following  Input 
(see  Appendix  A for  an  explanation  of  the  Input). 
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DYDIM  INPUT 


$ N30  - 61 , M30  ■ 61 , NP  - 1 $ 

NORMAL  INPUT  (8F10.2) 


0.14 

0.14 

0.18 

0.0 

2.8 

3.0 

0.6 

0.6E-06 

BLANK 

BLANK 

0.6 

100000. 

0.55 

5.55E-06 

24.0E-09 

0.16 

The  azimuthal  current  symmetry  around  the  cylindrical  load,  after  the  cur- 
rent flow  stabilizes.  Is  computed  to  be  29.3*.  The  current  flow  pattern  and 
azimuthal  current  symmetry  plots  for  this  run  are  shown  In  figure  11.  The  current 
symmetry  plot  Indicates  a current  asymmetry  of  about  30*  which  agrees  closely  with 
the  computed  value.  Unfortunately,  30*  asymnetry  may  not  produce  a viable  experi- 
ment so  an  attempt  must  be  made  to  Improve  the  current  symmetry. 

By  Increasing  the  Inductance  In  specific  regions  of  the  transmission  plate, 
the  current  symmetry  may  be  Improved.  One  such  example  Is  shown  In  figure  12.  This 
transmission  plate  Is  exactly  the  same  as  the  first  one  except  a wedge  on  each  side 
of  the  plate  has  higher  Inductance  by  a factor  of  10. 

This  higher  Inductance  may  be  obtained  by  milling  the  transmission  plates 
and  Inserting  additional  dielectric  material.  Thus,  the  capacitance  Is  also  de- 
creased by  a factor  of  10. 

The  Input  data  for  this  case  and 

OYDIM  INPUT 

$ N30  - 61 , M30  - 61 . NP  - 1 $ 

NORMAL  INPUT  (8F10.2) 

0.14  0.14  0.18  0.0  2.8  3.0 

0.6  0.6  100000.  5.55E-06  24.0E-09  0.16  .10  .10 

0.6E-06  0.55 

BLANK  CARO 

01213636 

01203737 

01193838 

01183939 

01174040 
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Figure  12.  Plate  Schematic  (With  Wedges). 


01164141 

01154242 

01144343 

01134444 

01124545 

01114646 

01104747 

01094848 

01084949 

01075050 

01065151 

01055252 

01045353 

01035454 

01025556 

01015760 

BLANK  CARO 

BLANK  CARO 

6Z7/8/9 


The  effect  of  the  wedges  Is  to  divert  the  current  so  that  It  flows  In  from 
the  four  comers  of  the  transmission  plate  rather  than  from  two  opposing  sides.  The 
current  flow  and  current  symmetry  plots  are  shown  In  figure  13.  (See  Appendix  F 
where  this  problem  was  used  for  sample  output.)  The  current  asymmetry  was  calcu- 
lated to  be  6.4%  at  0.18  m radius.  This  current  symmetry  should  be  adequate  for  a 
viable  experiment. 
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In  the  preceding  discussion,  two  geometries  were  compared  for  current 
symmetry.  The  second  system  was  Identical  to  the  first  except  for  the  Inclusion 
of  a high  Inductance  wedge.  How  much  does  the  wedge  cost  In  terms  of  an  Increase 
In  system  Inductance? 

PLATE  computes  partial  system  Inductance  at  specified  Intervals  of  time. 

This  Inductance  Is  reasonably  stable  as  a function  of  time  after  about  10  nsec. 
These  calculations  Indicate  an  Increase  In  total  system  Inductance  of  0.27  nH  due 
to  the  presence  of  the  wedges.  In  this  particular  example,  the  Inductance  In- 
crease Is  Intolerable  and  other  designs  must  be  tried.  One  possibility  Involves 
changing  the  Inductance  of  the  wedge  by  a factor  of  5 Instead  of  10.  A second 
possibility  Involves  decreasing  the  size  of  the  plate  but  leaving  the  wedge.  One 
such  design  will  provide  the  optimum  compromise  for  the  experiment  considered  here. 
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SECTION  V 

THE  CROSSED-PLATE  DESIGN 

Ourlng  the  development  of  a fast  one  megajoule  capacitor  bank,  the  crossed- 
plate  transmission  line  was  proposed.  A schematic  of  this  design  Is  shown  In  fig- 
ure 14.  Twenty  capacitor  modules  are  connected  to  the  four  "arms"  of  the  trans- 
mission plate. 

The  current  symmetry  of  such  a system  was  In  question;  consequently,  the 
capability  for  computing  such  a design  was  Incorporated  Into  PLATE.  This  option 
Is  called  by  setting  NSETUP  ■ 5 and  Introducing  various  Input  parameters  for 
capacitor  module  placement.  Typical  Input  for  the  problem  discussed  In  this  sec- 
tion Is: 

DYDIM  INPUT 

$ N30  « 780  , M30  - 109  , NP  ■ 1 $ 

NORMAL  INPUT  Format  (8F10.2) 


0.10 

0.14 

0.20 

0.0 

2.8 

5.0 

1.0 

7.78 

1.07 

100000. 

5.55E-06 

24.0E-09 

1.0 

13.33 

0.  IE-06 

10000. 

1.75 

2.542 

4.542 

5.152 

7.152 

BLANK 

BLANK 

This  calculation  yielded  the  current  path  plots  shown  In  figure  15. 

The  second  plot  Is  a magnified  view  of  the  first,  looking  only  In  the  vi- 
cinity of  the  load.  (There  Is  no  relative  current  density  associated  with  the 
current  flow  lines.) 

The  current  asymmetry  was  computed  to  be  1.5X.  Figure  16  shows  the  current 
as  a function  of  azimuthal  angle  around  the  load.  The  first  plot  Is  unsmoothed 
and  shows  the  effect  of  granlness  dur  to  the  cell  size  (In  this  case  the  cells  are 
1 cm  square).  The  second  plot  Is  smoothed  over  a 10  degree  azimuthal  angle,  and 
It  agrees  well  with  the  calculated  asymmetry.  This  current  symmetry  Is  definitely 
sufficient  for  a successful  experiment. 

The  load  Is  approximated  Inductively;  however,  Instead  of  resistive  losses, 
the  current  Is  deposited  In  large  capacitors  after  passing  through  the  simulated 
load.  One  drawback  Is  that  the  load  Is  considered  to  have  a static  rather  than 
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variable  Inductance.  The  capacitor  banks  and  transmission  plate  are  all  given 
realistic  values  of  Inductance  and  capacitance. 

The  above  considerations  Indicate  that  one  should  get  a realistic  time  de- 
pendent current  trace  and  calculated  one-eighth  system  Inductance.  Figure  17 


Figure  17.  Computed  Current  Trace. 


shows  a 2.32  usee  smoothed  current  trace  calculated  with  a coarser  mesh  than  was 
previously  considered.  The  current  peak  Is  at  1.16  usee  as  called  for  In  the 
specifications.  The  maximum  bank  current  Is  just  eight  times  the  current  In  the 
graph,  or  17.12  MA.  The  specifications  call  for  a maximum  current  of  about  20  MA. 

The  calculated  Inductance  per  one-eighth  system  Is  42.7  nH.  This  value  Is 
divided  by  eight  to  yield  the  total  system  Inductance  of  5.34  nH.  The  specifica- 
tions call  for  a system  Inductance  of  2 nH  outside  a meter  radius.  The  above 
PLATE  calculation  was  done  with  a load  at  a radius  of  20  cm.  The  formula 

L ■ UgAZ  in  (Ra/Ri) 
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may  be  used  to  compute  the  difference  In  Inductance  between  20  cm  and  100  cm  radii. 
(R2  ♦ 100  cm,  Rj  ■ 20  cm,  and  &z  ■ 0.0015  m.)  The  difference  Inductance  Is  3.03 
nN.  The  proposed  system  Inductance  for  this  load  Is  thus  5.03  nH.  The  calculated 
Inductance,  again.  Is  5.34  nH. 

The  two  Inductances  agree  very  well,  considering  the  approximations  made. 
(The  Inherent  Inaccuracy  In  the  Inductance  calculation  Is  dlsucssed  In  Appendix  0.) 
Further  error  arises  from  the  coarse  mesh  needed  to  run  the  crossed  transmission 
plate  problem  for  2 usee  of  real  time.  The  coarse  mesh  grossly  approximates  the 
load. 

The  system  Inductance,  current  profile,  and  maximum  current  are  close  to 
the  system  specifications.  The  design  of  this  1 mega joule  system  Is  thus  vali- 
dated by  PLATE. 
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SECTION  VI 
CONCLUSIONS 

1.  PLATE  uses  a simple  numerical  scheme  In  two  dimensions  to  simulate  the 
current  flow  In  parallel  plate  transmission  lines. 

2.  PLATE  Is  limited  to  five  transmission  plate  designs.  One  of  these  In- 
volves a one-dimensional  transmission  line.  The  other  four  are  two-dimensional 
problems  with  either  four  or  eight  fold  symmetry. 

3.  When  proper  zoning  Is  used  (usually  as  fine  as  1 centimeter  square 
zones),  the  measured  asymmetry  around  the  circumference  of  a cylindrical  load  Is 
about  It  accurate  at  a radius  of  30  centimeters. 

4.  PLATE  Is  a useful  tool  for  designing  transmission  plates  where  current 
symnetry  to  a cylindrical  load  Is  a valid  consideration. 

5.  Current  syimnetry  around  a cylindrical  load  may  be  Improved  by  current 
"shaping"  techniques.  One  such  technique  Is  the  Inclusion  of  high  Inductance 
areas  In  the  transmission  plates. 

6.  The  one  megajoule  crossed  plate  transmission  line  will  yield  current 
asynmetrles  less  than  2t.  The  system  Inductance  characteristics  are  verified  by 
PLATE  within  the  errors  Imposed  on  this  calculation. 
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APPENDIX  A 

DEFINITIONS  OF  INPUT  VARIABLES 


CARO  1: 
R1 : 

R2: 

RAD: 


R3: 

ER: 


SETUP: 


SETUP2: 


SETUP3: 


FORMAT  (8F10.2) 

Radius  of  load 

Radius  of  Increased  plate  separation  (milling) 

Radius  or  position  of  current  assymmetry  calculation  for 
plotting  and  first  of  15  radii  for  calculating  current 
assymmetry 

Radius  of  capacitor  bank  for  synmetry  test  problem 

Relative  permittivity  of  Insulator  between  transmission 
plates 

(1-5)  Determine  which  type  of  problem  will  be  run. 

1 Symmetry  test  problem 

2 Transmission  Line  Problem  with  matching  side 
boundary 

3 Rectangular  Transmission  Plate  with  two  capacitor 
modules 

4 Square  Transmission  Plate  with  four  capacitor 

modules 

5 Crossed  plate  transmission  line 

If  SETUP2  f 0,  the  program  Is  set  to  run  for  a long  time  and 
to  pick  off  I vs  time  plots  at  positions  halfway  between 
DIST1  and  DIST2  and  halfway  between  DIST3  and  DIST4.  This 
also  makes  a current  trace.  Use  only  when  SETUP  ■ 5. 

* 0,  nothing,  i 0,  the  value  of  cell  Inductance  In  the  simu- 
lated load. 


CARD  2:  FORMAT  (8F10.2) 

X:  X dimension  of  transmission  plate 

Y:  Y dimension  of  transmission  pltte 

Vo:  Voltage  on  capacitor 

CT:  Capacitance  of  1/2  module 

ALT:  Inductance  of  1/2  module 

CLENGTH:  Physical  length  of  1/2  capacitor  bank  module 

FACTOR:  Ratio  of  annular  plate  separation  to  normal  plate  separation 

FACT0R2:  Not  used  until  read  In  on  Card  5 

. 
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CARO  3:  FORMAT  (8F10.2) 

| TSTOP:  Time  In  sec  when  program  terminates 

DELZZ:  Plate  separation  (Default  ■ .0015m) 

OISTO:  X position  In  plate  where  current  Is  summed  (X  < DIST0 

< DELX,  or  error  mode  1.  results) 

DIST1:  First  of  four  variables  used  to  position  capacitor  banks 

0IST2:  All  four  are  the  boundaries  of  the  modules  on  the  side 

DIST3:  Of  the  crossed  plate  geometry.  In  the  case  of  the  trans- 

mission 

0IST4:  Line  OIST  1 and  DIST  2 position  the  capacitor  bank 

CARO  4:  FORMAT  (8F10.2) 

TEDIT  or  A:  Initial  edit  time  In  sec.  If  zero  default  is  TEOIT  » 40  nsec 

OELEO  or  B:  Increment  for  additional  edits.  If  zero  default  Is  DELED  - 

20  nsec 

TEDIT2  or  0:  Secondary  edit  time  and  Increment  In  sec,  If  zero  default  Is 

TEDIT2  - 10  nsec 

INPUT  TO  PLATE 

All  real  Input  Is  In  the  format  8F10.2 
All  Integer  Input  Is  In  the  form  512 


R1 

R2 

RAD 

R3 

ER 

SETUP 

SETUP  2 

SETUP  3 

X 

Y 

V 

CT 

ALT 

CLENGTH 

FACTOR 

FACTOR  2 

TSTOP 

TEDIT 

DELZZ 

DELED 

DISTO 

TEDIT2 

DIST1 

DIST2 

DIST3 

DIST4 

M N K L NQ 

FACTOR  2 (Only  If  NQ  f 0) 

M N K L NQ 

• 

• 

M N K L NQ 

FACTOR  2 (Only  If  NQ  f 0) 

BLANK  (To  terminate  Input) 

BLANK 

Input  terminates  when  K ■ 0 

NOTE:  FACTOR  2 Is  read  In  after  M N K L NQ  only  when  NQ  f 0.  Otherwise,  the 
previous  value  of  FACTOR  2 Is  used  for  Inductance/capacitance  changes. 
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NOTE  2:  H,  N,  K,  and  L are  limits  on  the  I and  J subscripts  which  define  the  X 
and  Y boundaries  within  which  the  Inductance  and  capacitance  are  changed 
by  the  factor  ■ FACT0R2.  K and  L are  the  lower  and  upper  limits  on  I, 
respectively,  and  they  thus  define  X boundaries.  M and  N are  the  lower 
and  upper  limits  on  J,  respectively,  and  they  thus  define  Y boundaries. 

The  only  additional  tnput  Involves  the  variable  dimensioning  feature  Included  by 
DYDIM.  This  Input  Is  of  the  form: 

$ N30  - 145,  M30  • 22,  NP  ■ 1 $ 

This  card  Is  usually  Included  or  changed  using  update.  N30  and  M30  are  the 
X and  Y dimensioning  of  the  arrays.  NP  Is  the  array  size  for  time  dependent  cur- 
rent plots  selected  when  SETUP2  f 0.  If  SETUP2  ■ 0 select  NP  ■ 1. 

CAUTION:  N30  and  M30  should  be  selected  so  that  DELX  and  DELY  are  equal;  other- 
wise, current  will  not  flow  properly  In  the  simulated  transmission  plate. 
DELX  and  DELY  are  determined  from: 

DELX  - X/(N30-2)  and  DELY  - Y/ (M30-2) 
where  X and  Y are  the  plate  dimensions. 


APPENDIX  B 
LIST  OF  VARIABLES 


A list  of  Important  variables  Is  found  below.  In  some  cases  the  variables 
have  a real  and  Integer  representation  because  they  are  read  In  as  real  variables 
and  used  as  Integers. 

A 

Initial  edit  time  ■ TEDIT 

AINCX 

Number  of  Increments  In  the  X direction 

AINCY 

Number  of  Increments  In  the  Y direction 
AJH 

Horizontal  current  array 

AJHQ 

ECS  array  of  AJH 
AJV 

Vertical  current. array 

AJVQ 

ECS  array  of  AJV 
AK 

No  longer  used 
AL 

Inductance  array 
ALN 

No  longer  used 
ALQ 

No  longer  used 

ALQQ 

ECS  array  of  AL 
ALT 

Total  Inductance  of  1/2  module 
AL1 

Inductance  of  each  simulated  capacitor  In  the  bank 
AL2 

Inductance  of  a normal  cell 

AMUO 

Permeability  of  free  space 
ANU 

No  longer  used 
AT 

Time  array  for  plotting 
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B 

Time  Increment  for  addition  edits  • DELED 
C 

Capacitance  array 
CLENGTH 

Length  of  1/2  module  of  capacitor  bank 
CQ 

ECS  array  of  C 
CT 

Total  capacitance  for  1/2  module 

CTRAC 

Array  for  plotting  time  dependent  current  trace 

CTR1 

Array  for  plotting  I trace  of  one  side  module  (one  megajoule  bank) 

CTR2 

Array  for  plotting  I trace  for  other  side  module  (one  megajoule  bank) 

Cl 

Capacitance  of  one  simulated  capacitor  In  bank 
C2 

Capacitance  of  a normal  cell 
C3 

No  longer  used 
C4 

Capacitance  of  a load  cell 
D 

Secondary  edit  time  • TEDIT2 

DELED 

Edit  Interval 

DELT 

Time  Increment 

DELX 

X dimension  of  cells 

DELY 

Y dimension  of  cells 

OELZ 

Normal  plate  separation 

OELZZ 

Input  variable  for  OELZ 

OISTO 

X position  In  plate  where  current  Is  summed 
DIST1  - DIST4 

The  four  boundaries  of  the  two  capacitor  modules  on  the  sides  of  the  one 
megajoule  transmission  line  arms. 

EPSI 

c.  permittivity  of  the  dielectric 
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EPS  10 

eQ,  permittivity  of  free  space 
ER 

Relative  permittivity  of  the  dielectric 
FACTOR 

Ratio  of  annular  plate  separation  to  normal  plate  separation 
FACT0R2 

Ratio  of  special  regions  to  normal  plate  separation 
I 

Do  loop  Index  In  X direction 
III 

No  longer  used 

INCX 

Same  as  AINCX 

INCXO 

INCX  - 1 

INCX1 

INCX  ♦ 1 

INCY 

Same  as  AINCY 

INCY1 

INCY  + 1 
K 

Index  for  the  major  loop  over  time  also  passed  from  WEDGE  to  SETUPX  as  a 
plotting  array  size 

KSEN 

No  longer  used 
KSENSW 

No  longer  used 

NCAP 

Number  of  capacitors  In  one-half  module 
NDIST1  - NDIST4 

Corresponds  to  DIST1  - DIST4  except  these  are  cell  designations  of  the 
boundaries  of  the  two  capacitor  modules 

NN 

NDISTl  or  NDIST3 
NNN 

NDIST2  or  NDIST4 

NECS 

Decimal  ECS  needed 
NP0SX1 

Position  of  one  I trace 
NP0SX2 

Position  of  second  I trace 
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NSETUP 

Type  of  problem  being  run 
NSETUP2 

Cross  plate  only  - picks  off  I traces  as  well  as  I vs  t and  plots  them 
NSETUP3 

If  I1  0 read  In  as  the  value  of  cell  Inductance  In  the  simulated  load 
NT 

Number  of  time  Increments 

NZAP 

An  Integer  variable  used  for  filling  plotting  arrays 

N100 

Subscript  of  plotting  arrays  while  they  are  being  filled. 

RAO 

Radius  where  current  assymnetry  Information  Is  tabulated 

RES  IS 

Resistance  of  one  cell 
R1 

Radius  of  load 
R2 

Outside  radius  of  annular  region  of  Increased  plate  separation 
R3 

Radius  of  Inner  edge  of  capacitor  bank  when  symmetry  test  problem  Is  run 
R4 

No  longer  used 
T 

Problem  time 

TEDIT 

Initial  edit  time 
TEDIT2 

Secondary  edit  time  and  Increment 

TSTOP 

Stop  time  of  problem 

V 

Voltage  array 
VO 

Initial  capacitor  bank  voltage 
VQ 

ECS  array  of  V 

VI 

Scratch  array 
HA 

Scratch  array 
WB 

Scratch  array 
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APPENDIX  C 

MUTUAL  INDUCTANCE  EFFECTS 


The  program  PLATE  Ignores  mutual  Inductance.  This  appendix  justifies  that 
assumption  and  gives  a feeling  for  the  magnitude  of  error  Introduced.  The  mutual 
Inductance  between  adjacent  cells  Is  calculated,  assuming  the  cells  are  like 
small  Inductors  In  parallel,  with  the  accompanying  fringing  which  leads  to  this 
mutual  Inductance.  It  Is  concluded  that  mutual  Inductance  has  only  a small  effect 
on  total  plate  Inductance  as  long  as  the  cell  size  Is  large  compared  to  the  plate 
separation. 

Figure  C-1  shows  two  adjacent  cells  where  a discontinuity  In  plate  separa- 
tion exists  at  the  Interface.  Such  discontinuities  can  be  used  for  selectively 
controlling  cell  Inductance. 


Figure  C-l.  Adjacent  Cells. 

If  each  cell  Is  considered  to  be  a current  loop,  the  mutual  Inductance  of  a 
half-loop  acting  on  a whole  current  loop  as  shown  In  figure  C-2  can  be  determined. 

The  currents  In  the  Z direction  are  cancelled  by  adjacent  current  loops  and 
can  thus  be  Ignored.  The  flux  created  by  wire  A which  passes  through  current  loop 
B Is  given  by 


ylf  • n da 


(C-l) 
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u0i.t  r.  ^,vr  r.  i^vr'1 

*ab " u [1+  v — sw  / J [1+  {— ary  J (c-4) 

For  the  bottom  half  of  current  loop  A,  la  Is  negative  and  the  limits  on  r 
are  reversed;  hence,  there  are  two  equal  contributions  to  the  flux.  Since 

AD 


ab/Ia’ 


1/2  r M-^vr1'2 


y.i  r /az2+azA2"]  r / 

■ -r  ^ ['*  Htr-j  J H 


In  general  the  self-inductance  Is  given  by 


L " uo  IT 


and  the  coefficient  of  coupling  Is  given  by 


k‘r 


Substituting  for  M and  L the  coefficient  of  coupling  becomes 


For  the  specific  case  where  aZ2  ■ 1lAZlt  and  W ■ lOAZi, 


k - li  tnJl+JJ/*)2  . 1°  (0.04217)  - 0.134 
w v l+(l/2)2  w 


(C-7) 


(C-8) 


(C-9) 


Thus  the  mutual  Inductance  Is  13.4*  of  tfte  smaller  Inductor  and  about  1.2*  of  the 
larger  Inductor. 

For  cases  where  aZj  ■ aZ2,  the  coefficient  of  coupling  reduces  to 


k * *a7 
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; 


I 


when  aZ  < W,  the  In  term  nay  be  expanded  so  that 

k " *aT  (1/2)  [(t)2  ’ (1/2)  + C1/3)  H)S  ‘ * * *] 

or 

k’tM’  • (V2)  * (l/3)  H)'  * ' • •] 

wheyi  aZ  « U only  the  first  term  Is  needed  so  that 


k 


aZ 

lirU 


(C-12) 


This  coefficient  of  coupling  between  cells  with  plate  separation  aZ  - 0.1 
AX  Is  0.016.  Clearly,  all  mutual  Inductance  effects  can  be  kept  Insignificant  by 
keeping  the  ratio  of  aZ  to  aX  smaller  than  0.1.  This  limit  Is  approximately  main 
talned  In  PLATE  calculations  except  In  the  simulated  load  and  the  high  Inductance 
wedges,  where  the  plate  separation  and  cell  size  are  approximately  equal. 

A violation  of  the  above  limit,  where  high  Inductance  wedges  are  Included, 
Is  discussed  In  section  IV.  For  that  calculation,  the  plate  separation  Is  1.5 
times  the  cell  size  In  a wedge  shaped  area.  The  total  Increase  In  system  Induc- 
tance Is  only  0.27  nH.  The  mutual  Inductance  effect  on  the  system  Inductance  Is 
maximized  near  this  discontinuity.  The  coefficient  of  coupling  of  a cell  at  the 
discontinuity  Is  calculated  [using  equation  (C-10)]  to  be  0.58 7 when  compared  to 
a normal  cell  and  0.0587  when  compared  to  a wedge  cell.  Fortunately,  the  wedge 
Interface  Into  the  plate  Is  over  a limited  region.  The  effective  mutual  Induc- 
tance of  the  wedge  Is  expected  to  add  less  than  1.0  nH  to  the  system. 

Thus  the  typical  transmission  plate  Inductance  Is  accurate  to  within  5%  If 
the  criterion  aZ*£0.1aX  Is  met  everywhere  except  In  the  load  and  the  high  Induc- 
tance wedges.  If  this  criterion  Is  met  everywhere,  the  transmission  plate  Induc- 
tance Is  accurate  to  within  1.6X. 


AFWL-TR- 76-205 


APPENDIX  D 
VALIDATION  OF  PLATE 


Two  test  problems  are  solved  In  this  Appendix.  The  first  Is  a transmis- 
sion plate  with  Imposed  azimuthal  symmetry.  The  calculations  show  that  current 
Is  propagated  with  azimuthal  symmetry  through  the  square  mesh  and  that  current 
symmetry  at  the  simulated  load  Is  excellent.  The  second  problem  Is  a parallel 
plate  transmission  line  with  a capacitor  bank  on  one  side  and  a load  on  the 
other.  The  system  Inductance  Is  calculated  two  ways:  In  PLATE  and  analytically 
using  equation  (7)  In  section  IV.  The  two  results  agree  nicely.  The  PLATE  calcu- 
lation was  accomplished  with  two  different  cell  sizes.  Cell  size  Is  shown  to 
have  only  a small  effect  on  system  Inductance. 

1.  Symmetry  Test  Problem 

This  problem  Is  a specific  option  In  PLATE.  It  Is  called  with  SETUP  ■ 1. 
The  Input  cards  for  this  problem  with  25  x 25  zones  are 

DYDIM  Input 

$ N30  - 27,  M30  - 27,  NP  • 1 $ 

Normal  Input 


0.259 

0.406 

0.3 

0.8 

2.8 

1.0 

1.0 

0.6E-07 

BLANK 

• 

BLANK 

1.0 

100000. 

1.0 

5.55E-06 

24.0E-09 

1.0 

10.0 

The  current  flow  pattern  and  current  symmetry  plots  are  shown  In  figure  0-1. 
The  minimum  calculated  current  asymmetry  was  5.6X  (the  current  asymnetry  Is  calcu- 
lated at  many  radii  outside  the  load).  The  cells  had  a size  of  4 cm  square.  This 
coarse  mesh  contributes  to  the  current  asymmetry. 

A second  problem  was  run  with  100  x 100  zones.  The  cell  size  was  1 cm 
square.  The  minimum  calculated  current  asymmetry  was  IX.  This  IX  current  asym- 
metry Is  excellent,  considering  It  was  calculated  around  a quarter  circle  that  was 
superimposed  on  square  cells.  The  current  flow  equations  are  shown  to  be  spatially 
valid,  since  current  does  not  flow  preferential ly  In  the  horizontal,  vertical,  or 
diagonal  direction  relative  to  the  square  grid. 
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Th«  current  flow  pattern  and  current  symmetry  plots  for  this  run  are  shown 
In  figure  0-2. 

These  two  plots  agree  In  substance,  but  the  first  shows  the  effect  (graini- 
ness) of  the  relatively  large  4 cm  square  cells.  Current  flow  varies  drastically 
with  angle  at  the  plotted  radius.  The  second  plot  shows  an  almost  constant  cur- 
rent flow  for  different  angles. 

2.  Parallel  Plate  Transmission  Line 

The  second  test  problem  Involves  a parallel  plate  transmission  line  with 
matching  side  boundaries.  This  approximates  a cylindrical  co-axial  transmission 
line. 


A transmission  line  of  this  sort  may  be  represented  as  an  Inductor  con- 
nected between  two  capacitors.  One  of  the  capacitors  Is  the  capacitor  bank,  and 
It  Is  charged  to  a predetermined  voltage.  The  second  capacitor  Is  a short  cir- 
cuit load. 


The  current  In  the  Inductor  Is  given  by 


41-  Y 

ar  r 


or 


(D-l) 


t* 


At  AV 
^n+l  " 1 


(0-2) 


where  V Is  the  average  voltage  across  the  Inductor,  t Is  the  time  Interval, 

(Jn+1  _ V 1*  the  Increase  In  the  current  during  this  time  Interval,  and  L Is  the 
effective  Inductance  of  the  system.  This  equation  Is  used  to  determine  the  Induc- 
tance of  the  system. 


The  accuracy  of  this  method  Is  limited.  There  Is  a wave  superimposed  on  the 
plotted  current  trace.  This  wavy  nature  Is  caused  by  waves  sloshing  around  the 
simulated  plate  because  the  capacitor  bank  and  load  are  not  exactly  matched.  Be- 
cause of  these  waves,  relatively  long  Intervals  of  time  must  be  considered  to  make 
the  Inductance  calculations  meaningful. 
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The 

transmission  line  problem  Is  specified  with 

the  following  Input: 

0.0 

0.0 

0.0 

0.0  2.8 

2.0 

1.0 

1.5E-07 

BLANK 

BLANK 

1.0 

0.0015 

100000. 

0.05 

5.55E-06  4.0E-09 

1.0 

OYDIM  Input  Is  also  used  to  specify  25  x 25  zones  or  50  x 50  zones.  The  Impor- 
tant part  of  the  above  Input  Is  the  capacitor  bank  Inductance  L ■ 4 nH. 

The  25  x 25  zones  and  50  x 50  zones  problems  yielded  system  Inductances 
6.15  nH  and  6.05  nH,  respectively,  between  50  nsec  and  100  nsec  problem  time 
after  Initiation. 

The  bank  has  an  Inductance  of  4 nH.  There  Is,  thus,  approximately  2.1  nH 
In  the  transmission  line. 

Using  equation  (9)  from  section  III,  the  Inductance  of  a parallel  plate 
transmission  line  with  dielectric  thickness  0.0015  m,  length  x ■ 1.0  m,  and  width 
y ■ 1.0  m Is  given  by 

L " “o  aT  42  " 4*  x 10"7  °*0015  * ]-88  nH 

This  value  Is  about  10*  less  than  the  value  of  Inductance  calculated  by  PLATE,  even 
so  this  comparison  helps  to  validate  the  numerical  method  used. 

The  plotted  current  trace  Is  an  approximation,  since  the  Inductance  calcu- 
lations art  tied  directly  to  the  current  values.  In  Z-pInch  type  problems,  the 
Inductance  of  the  load  Increases  as  the  foil  Implodes.  No  attempt  has  been  made  to 
simulate  the  rising  Inductance  In  these  calculations.  Thus,  the  Inductance  calcu- 
lations and  current  trace  do  not  approximate  the  real  experiments  after  0.3  usee 
(that  Is  after  the  foil  begins  to  Implode). 
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APPENDIX  E 
LISTING  OF  PLATE 
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PROGRAM  PLATE  < INPUT  .OUTPUT .MF35PL)  PLATE 

PLATE 

PLATE 

PLATE  IS  A COOE  WHICH  OETERMJnES  ThE  PATHS  ALONG  which  PLATE 

CURRENTS  FLO 'W  THROUGH  TRANSMISSION  PLATES.  A PLATE  IS  DIVIDED  PLATE 
INTO  A SQUARE  MESH.  EACH  CORNER  PLATE 

IN  The  mesh  HAS  A CAPACITOR  CONNECTED  TO  GROUNO.  The  VOLTAGE  FOR  PLATE 
The  next  TIME  STEP  ON  THIS  CAPACITOR  IS  DETERMINES  BY  THE  PLATE 

CURRENTS  IN  THE  FOUR  LINES  COMING  TO  THAT  CAPACITOR  BY  PLATE 

PLATE 

V2  » VI  - OELT»  < J1 ♦ J2» J3» J4) /C  PLATE 

PLATE 

WHERE  V2  IS  ThE  NEW  VOLTAGE  PLATE 

VI  IS  THE  OLD  V0LTA6E  PLATE 

OELT  IS  The  TIME  STEP  PLATE 

Jl»  J2*  J3.  J4  ARE  THE  FOUR  CURRENTS  WHICH  CONTAIN  PLATE 

INTERNAL  SIGNS  AND  wHERE  The  SIGN  CONVENTION  IS  IMPORTANT  PLATE 
C IS  THE  CAPACITANCE  PLATE 

PLATE 

The  links  IN  The  MESH  WHICH  CONNECT  The  CAPACITORS  CONTAIN  PLATE 
INOUCTORS  ANO  RESISTORS  IN  SERIES.  HE  CURRENTS  which  PROPAGATE  PLATE 
THROUGH  THESE  LINKS  ARE  CHANGEO  IN  THE  BY  PLATE 

PLATE 

J2  ■ IVA  - VB  * Jl* (L/OELT-R/2) ) / (L/OELT.R/2)  PLATE 

PLATE 

WHERE  J2  IS  THE  NEW  CURRENT  PLATE 

Jl  IS  ThE  OLD  CURRENT  PLATE 

VA.VB  ARE  ThE  VOLTAGES  ON  ThE  CAPACITORS  AT  EITHER  END  PLATE 
Of  The  INDUCTOR  ANO  RESISTOR  PAIRS.  AGAIN  The  SIGN  PLATE 

CONVENTION  IS  IMPORTANT.  WITH  SIGNS  CONTAINED  INTERNALLY  PLATE 
IN  VA  ANO  Vb.  PLATE 

L IS  The  INOUCTANCE  IN  A CONNECTING  LINK  PLATE 

R IS  The  RESISTANCE  IN  A CONNECTING  LINK  PLATE 

PLATE 

THE  main  ROUTINE  PLATE  IS  THE  CONTROLLING  ROUTINE  WHICH  PLATE 

CALLS  ThE  VARIOUS  SUBROUTINES.  SUBROUTINES  C8ANK  ANO  CriCw  PLATE 

CONTROL  CURRENT  FLOW  OUT  Or  THE  CAPACITOR  BANKS  ANO  THROUGH  THE  PLATE 
PLATE.  RESPECTIVELY.  OuRING  EACH  THE  STEP  The  CURRENT  IS  HOVEO  PLATE 
THROUGH  The  ENTIRE  mesh  USING  OLO  VOLTAGES.  Then  HE  VOLTAGES  PLATE 
ARE  CHANGED  USING  THE  NEW  CURRENTS.  PLATE 

THE  PROGRAM  BEGINS  AT  T«0  WITH  HE  CAPACITOR  BANK  BEING  PLATE 

SWITCFEO  ON.  ONE  MUST  NOTE  rt£«£  THAT  USUALLY  A SYMMETRIC  ONE-  PLATE 
QUARTER  OR  ONE-EIGHTH  OF  HE  PROBLEM  IS  CONSIDERED  OuRING  A RUN.  PLATE 
This  ALLOWS  A finer  HCSh  ANO  A MORE  ACCURATE  SOLUTION.  PLATE 

HE  PROGRAM  PLATE  CALLS  ThE  FOLLOWING  SUBROUTINES  PLATE 

THESE  ROUTINES  ARE  LISTEO  IN  THE  OROER  They  APPEAK  IN  The  PLATE 

listing.  PLATE 

PLATE 

CBANK  - PROPAGATES  CURRENT  OUT  OF  ThE  CAPACITOR  BANKS  ANO  PLATE 
INTO  THE  PLATE.  PLATE 

PLATE 

CFLOw  - PROPAGATES  CURRENT  THROUGH  THE  MESH  OF  THE  PLATE.  PLATE 

PLATE 


1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

H 

IS 

IS 

17 

IS 

19 

20 
21 
22 
23 
2“ 

25 

26 
27 
2S 

29 

30 

31 

32 

33 
3“ 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

4? 

46 

49 

so 

51 

52 

53 

54 
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MISCELL  - THIS  SUBROUTINE  READS  THE  INPUT  AND  PRINTS  IT  OUT. 
IT  uses  THIS  INPUT  TO  SET  UP  TH£  PROBLEM.  IT  DETERMINES  ThE  TIME 
STEP  ANO  OTHER  VARIABLES  NEEOEO  TO  RUN  The  PROBLEM.  Th£n  IT 
CALLS  ThE  FOLLOWING  THREE  SUBROUTINES  WHICH  ACTUALLY 
INITIATE  The  ARRAYS  NEEUEO  TO  RUN  THE  PROBLEMS. 

MESH  - THIS  SUBROUTINE  INITIALLIZES  THE  INDUCTANCE. 
CAPACITANCE.  VOLTAGE.  ANO  CURRENT  ARRAYS  TO  GENERAL  PLATE  VALUES. 

LOAO  - THIS  SUBROUTINE  RESETS  ARRAY  VALUES  TO  SET  UP  A 
SIMULATED  LOAD. 

WEDGE  - ThIS  SUBROUTINE  RESETS  ARRAY  VALUES  IN  INPUTEO 
REGIONS  OF  The  plate  TO  nCLR  IN  OBTAINING  BETTER  CURRENT 
SYMMETRY. 

SETUPX  - PLOTS  ON  MICROFILM  The  PROBLEM  SETUP  BY  DRAWING 
LINES  WHICH  SEPARATE  DIFFERENT  REGIONS  IN  The  PLATE. 

PRINT1 (NEOIT)  - PRINTS  EDITED  VERSIONS  OF  ThE  VOLTAGE. 
CAPACITANCE.  INDUCTANCE.  VERTICAL  CURRENT  OR  HORIZONTAL  current 
ARRAYS.  WHICH  ARRAY  IS  PRINTE3  OEPENOS  ON  ThE  CALLING  VARIABLE 
NEOIT. 

PRINT2  - PRINTS  SMALLER  EDITED  ARRAYS  OF  VOLTAGE. 

HORIZONTAL.  ANO  VERTICAL  CURRENT. 

INOUCT  - THIS  subroutine  cohputcs  ANO  prints  The  current 
flowing  IN  THE  PLATE  Each  nano-second,  it  prints  charge 
conservation,  it  also  saves  the  current  at  eoual  time  intervals 
SO  THAT  a TI-E  OEPENOENT  CURRENT  trace  CAN  BE  MAOE  BY  subroutine 
CTRAC. 

INOUCT 2 - This  SUBROUTINE  computes  ANO  PRINTS  SMOOTHED 
SYSTEM  INDUCTANCE  FOR  SPECIFIC  TIME  INTERVALS. 

CTRACE  - This  SUBROUTINE  PLOTS  A time  DEPENDENT  CURRENT 
TRACE.  TjiIS  SUBROUTINE  ANO  THE  NEXT  ARE  NOT  CALLED  EXCEPT  «HEN 
SCTUP2  IS  NOT  EQUAL  TO  0 IN  WHICH  CASE  THE  PROBLEM  IS  SET  UP  TO 
R'JN  with  A GROSS  MESH  ANO  FOR  A LUNG  TIME. 

CTRAC2  - THIS  SUBROUTINE  COMPUTES  ANO  PLOTS  DI/DT  TERMS  FOR 
THE  capacitor  mooules. 

CURRENT  - COMPUTES  ANO  PLOTS  SEVERAL  CURRENT  PATHS  FROM  Th£ 
CAPACITOR  BANK  TO  THE  LOAD. 

UNWINO  - COMPUTES  ANO  PLOTS  THE  THETA  CURRENT  SYMMETRY  OF 
THE  CURRENT  GOING  TO  ThE  LOAO. 

linplot  - plots  the  oata  on  microfilm. 

LIST  OF  VARIABLES  IN  DIFFERENCE  EQUATIONS 


PLATE  55 
PLATE  56 
PLATE  57 
PLATE  SB 
PLATE  59 
PLATE  60 
PLATE  61 
PLATE  62 
PLATE  63 

plate  &a 

PLATE  65 
PLATE  66 
PLATE  67 
PLATE  6a 
PLATE  69 
PLATE  70 
PLATE  71 
PLATE  72 
PLATE  73 
PLATE  7<* 
PLATE 
PLATE 
PLATE  77 
PLATE  73 
PLATE  79 
PLATE  80 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE  87 
PLATE  83 
PLATE  89 
PLATE  90 
PLATE  91 
PLATE  92 
PLATE  93 
PLATE  9h 
PLATE  95 
PLATE  96 
PLATE  97 
PLATE  99 
PLATE  99 
PLATE100 
PLATElOi 
PLATEI02 
PLATEI03 
PLATEIOH 
PLATE105 
PLATElOb 
PLATEI07 
PLATE108 


75 

76 


61 

82 

83 

89 

85 

86 
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V - VOLTAGE 

AJV  - VERTICAL  CURRENT 

Ajh  - horizontal  CURRENT 
c - capacitance 

AL  - INOUCTANCE 
AL1  - SANK  INOUCTANCE 
RESIS  - RESISTANCE 
OELT  - TIME  STEP 

VI  - SCRATCH  arrat  useo  tor  voltage  or  current 

SIMILAR  VARIABLES  ENDING  IN  0 ARE  TWO  DIMENSIONAL  ECS  ARRAYS  AND 
THET  REPRESENT  The  SAME  quantities. 


PLATE109 

PLATEUO 

PLATEUI 

PLATE112 

PLATEUO 

PLATEIU 
PLATE1 15 
PLATE116 
PLATE1I7 
PLATE1I8 
PLATE119 
PLATE12C 

plate  121 
•plate i 22 
PLATE123 
PLATE12A 


COMMON  X. 

Rl. 

VO. 

INCX. 

NOISTl. 

OISTO. 

RESIS. 

AJV (26) « 

PLATE12S 

1 

Y. 

R2. 

CT. 

INCY. 

NO  I S T 2 . 

INCXO. 

*SENSM . 

AJH (26) . 

PLATE126 

2 

K. 

R3. 

ALT. 

AINCX. 

N0IST3. 

TSTOP. 

NSETUP. 

V (26) • 

PLATE127 

0 

T. 

*4. 

NT. 

AINCY. 

NOISTa. 

TEOIT. 

NSETUP2. 

C (26) . 

PLATE123 

4 

AKf 

RAO. 

ALN. 

INCX 1 . 

3IST1. 

DELED. 

NSETUP3. 

V 1 (26)  . 

PLATE129 

5 

CELT. 

AL1. 

ALQ. 

INCY 1 . 

0IST2. 

TE0IT2. 

NSETUPA, 

AL (26) • 

PLATE130 

6 

KSEN. 

AL2. 

ANU. 

UELX. 

0IST3. 

NPOSX I . 

NSETUPS. 

*4(26.26)  • 

PLATElOl 

7 

NCAP, 

NIOC. 

NZAP. 

oely. 

OISTa. 

NP0SX2. 

Ill  (10). 

•8 (26.26) • 

PLATE  1 32 

8 

XI  ( 100) 

. Y 1 ( 100) . 

AT  < 1 ) * 1 

CTRAC  < l ) . 

CTRH1) 

. CTH2U) 

PLATEUO 

DIMENSION  AJVQ (26f 26) t AJHO(26.2b> . VQ(26.26).  C3I26.26). 
1 ALQG (26*26) 

LEVEL  3.  AJVQ.  AJHQf  VQt  CO*  ALSO 

COMMON  /ECS1/  AJVQ 
COMMON  /ECS2/  AJHO 
COMMON  /EC SO/  VQ 
COMMON  /ECS*/  CQ 
COMMON  /ECSS/  ALQG 


eix  ecs  riELo  length 

CALL  INITPLT 
CALL  MISCELL 
T*0.0 
N 100-0 
KSENSVaO 

IE  OELX  DCES  NOT  EQUAL  DELY  THE  CALCULATION  TERMINATES. 

IE  C <0ELX-0ELY>/0ELX.LT.1.0/AMAXUAINCA. AInCY) ) GO  TO  1 
PRINT  7 
PRINT  8 
STOP 

CONTINUE 


PLATE10A 
PLATE135 
PL ATE  1 06 
PLATE137 
PLATE138 
PLATE109 
PLATEUO 
PLATEUI 
PLATEU2 
PLATEUO 
PLATEUa 
PLATEU5 
PLATEUO 
PLATEU7 
PLATE  U8 
PLATEU9 
PLATE150 
PLATE151 
PLATE152 
PLATE1S0 
PLATE1SA 
PLATEISS 
PLATE156 
PLATEIS7 
PLATE1S8 
PL4TE1S9 
PLATEUO 
PLATEUI 
PLATE162 


uo  uuuu  «m  n#  m our- 


r< 


8 


THESE  LOOPS  PERFORM  THE  ITERATIONS  OVER  A TINE  STEP  DELT. 

00  6 K»1 »NT 
T«T »0ELT 
CALL  CSANK 
CALL  CFLOw 

THE  REST  OF  THIS  PROGRAM  CALLS  VARIOUS  CHECKING.  PRINTING. 
ANO  PLOTTING  SUBROUTINES. 

IF  (NSETUP2.E0.0)  GO  TO  2 
IF  (NZAP.EQ.O)  GO  TO  3 
GO  TO  A 
CONTINUE 

IF  CT.LT.TEDIT2)  GO  TO  A 

TE0IT2«TEDIT2*i.0E-09 

CONTINUE 

continue 

IF  (T.LT.TE0IT)  go  to  s 

tedit«teoit.oeleo 

KSENSwaO 

CALL  PRINT2 
CALL  CURRENT 
CALL  UNWIN0 
CONTINUE 

CALL  ATIME  (CP.PP.IO.TIMTGO) 

IF  (TIMTOO.GT. 30.0. AND. T.LE.TSTOP)  GO  TO  6 
CALL  CTRACE 

IF  (NSETUP2.GT .0)  CALL  CTRAC2 
STOP 

CONTINUE 


PLATEI63 
PLATE 16A 
PLATE16S 
PL ATE  166 
PLATE167 
PLATE  168 
PLATE169 
PLATE170 
PLATE171 
PLATE172 
PLATE173 
PLATE17A 
PLATE! 75 
PLATE176 
PLATE177 
PLATE178 
PLATE179 
PLATE1S0 
PLATE181 
PLATE182 
PLATE183 
PLATE  1 8A 
PLATE185 
PLATE186 
PL ATE  1 87 
PLATE  1 88 
PLATE189 
PLATE ! 90 
PLATE191 
PLATE192 
PLATE193 
PLATE19A 


FORMAT  <5X,72HTHIS  PROGRAM  TERMINATED  BECAUSE  OELA  010  NOT  EQUAL  CPLATE195 


1ELT  (APPROXIMATELY).)  - PLATE196 

format  (5X.SAHCORRECT  OTOI**  INPUT  N]0  OR  M30  TO  ADJUST  OELA  OR  0ELPLATE197 


IT.)  ~ PLATE  193 

ENO  - PLATE199 


SUBROUTINE  C8ANK 


DIMENSION  AJVO (26*26) * AJhQ(26<26)«  VO(26*26)«  C0(26.26>* 
I ALQu (26*26) 

LEVEL  3*  AJVQ.  AJHO,  VO*  CO*  ALOO 

COMMON  /CCSI/  AJVO 
COMMON  /ECS2/  AJnO 
COMMON  /ECS3/  VQ 
COMMON  /ECSA/'  CO 
COMMON  /ECS5/  ALOO 


IF  NSETUP  OOES  NOT  equal  1 the  LOOPS  TO  STATEMENT  6 SEND  CURRENT 
INTO  THE  PLATE  FROM  TnE  CAPACITOR  BANKS. 

IF  (NSETUP. EO.l)  GO  TO  5 

I«INCX1 

CO  l J*1.NCAP 

CALL  REAOEC  ( AJV ( I ) * AJVO ( I * J) « 1 ) 

CALL  REAOEC  (V(l) .VQ(I-l.J) *2) 

CALL  REAOEC  (C(I) .CO(I.J).l) 

AJV(I)bAJV(D*0ELT*(V(2)-V(I))/AL1 
V (2)  »V ( 2) -AJV ( I ) »OELT/C ( I ) 

CALL  WRI TEC  (V (2) , VO ( I . J) . 1 ) 

CALL  nRITEC  (AJV ( I) ,AJVQ(I,J> *1) 

CONTINUE 

IF  NSETUP«5  CURRENT  IS  FEO  FROM  The  CAPACITOR  BANKS  ON  The  SICE 
OF  The  CROSSED  TRANSMISSION  PLATE.  Th£  TkO  SIOE  CAPACITOR 
MODULES  HAVE  LIMITS  OF  NUIST1*  N0IST2.  NOIST3.  Anu  NOISU. 

IF  (NSETUP. NE. 5)  GO  TO  5 

NNaNOISTl 

NNN.NO I ST2 

J-INCTl 

CALL  REAOEC  (AJM( 1) «AJHQ( I* J) • INCX) 

CALL  REAOEC  (V(1)*VQ(1*J). INCX) 

CALL  REAOEC  (V!(i) * VQ ( 1 * J— 1 ) * INCX) 

CALL  REAOEC  (C(I)*CU(1*J>* INCX) 

NZAPmMOO(K.IOO) 

IF  (NZAP.NE.O.OR.NSETUP2.EO.O)  GO  TO  2 
N10®.N1«9»1 


PLATE200 
PLATE201 
PLATE202 
PLATE203 
PLATE204 
PLATE205 
PLATE206 
PLATE297 
PLATE208 
PLATE209 
PLATE210 
PLATE21 I 
PLATE212 
PLATE213 
PLATE21A 
PLATE21S 
PLATE216 
PLATE217 
PLATE2I8 
PLATE219 
PLATE220 
PLATE221 
PLATE222 
PLATE223 
PLATE224 
PLATE225 
PLATE226 
PLATE227 
PLATE228 
PLATE229 
PLATE230 
PLATE231 
PLATE232 
PLATE233 
PLATE23A 
PLATE235 
PLATE236 
PLATE237 
PLATE239 
PLATE239 
PLATE2A0 
PLATE241 
PLATE242 
PLATE2A3 
PLATE2*<* 
PLATE24S 
PLATE26* 
PLATE2A7 
PLATE2A8 
PLATE249 
PLATE250 
PLATE251 
PLATE2S2 
PLATE253 


- ■ 


IF~~r 


CTR1(N100>»<V«NPOSX1)-VI<NPOSX1>>»DELT/AL1 
CTR2 (N100) ■ (V (NPOSX2) ”Vl (NPOSA2) ) *DELT/ALl 

2 CONTINUE 

00  3 I »NN  t NNN 

AJH«I>«AJH(I)MV<H-VHII  >*0€LT/ALI 
VU)«vat-OELT»AJH<l)/C«I) 

3 CONTINUE 
NN-N0IST3 
NNMN0IST4 

00  4 I «NN  • NNN 

AJHU)«AJM<I>»<V<I)-Y1(I>  >«oelt/ali 
V(I)«VU>-0ELT*AJH<U/C<I) 

4 CONTINUE 

CALL  'API TEC  (VII)tVQ(l.J)tlNCX) 

CALL  WRITEC  < A JH i I ) • A JHO ( 1 • J ) • I NCX ) 

5 CONTINUE 
RETURN 
ENO 


PLATE2S4 

PLATE255 

PLATE256 

PLATE257 

PLATE253 

PLATE259 

PLATE260 

PLATE261 

PLATE262 

PLATE263 

PLATE264 

PLATE263 

PLATE266 

PLATE26T 

PLATE268 

PLATE269 

PLATE270 

PLATE271 


SUBROUTINE  CFLOW 


PLATE272 

PLATE273 


COMMON  X* 

Rl* 

vo* 

INCX. 

NOISTl. 

OISTO* 

RESIS, 

AJV (26) . 

PLATE27* 

1 

Y* 

R2. 

CT. 

INCY. 

N0IST2. 

INCXO. 

K SENS A. 

AJH (26) t 

PLATE275 

2 

K* 

R3. 

ALT* 

A INCX. 

MOIST  3* 

TSTOP. 

NSETUP, 

V (261 • 

PLATE276 

3 

T. 

P*. 

NT. 

AINCY. 

MOIST*. 

TEOIT. 

NSETUP2. 

C (26) • 

PLATE277 

* 

AK  • 

RAO* 

ALN* 

INCX 1 * 

OIST 1 * 

OELEO* 

NSET0P3. 

VI (26) . 

PLATE279 

S 

OELT. 

AL1. 

ALO* 

1 MCY 1 . 

01ST2, 

TE0IT2. 

NSETJP*. 

AL (26) * 

PLATE279 

6 

KSENt 

AL2* 

ANU* 

oelx  . 

0IST3, 

NP0SX1. 

NSETUP5. 

«A(26.26> . 

PLATE280 

7 

NCAPf 

NIOOt 

NZAP* 

UCLY. 

OIST*. 

NP0SX2. 

IIIUO). 

»o 126.26) . 

PLATE281 

8 

XI  ( 180)  * YldOO). 

ATI l)  * 

CTRAC(l) 

« CTRld) 

. CT*2 ( 1) 

PLATE282 

01 PENSION  AJVQ (26.26)  » AJHQI26.26).  VQI26.26).  CU<26.26>. 
1 ALOC (26*26) 

LEVEL  3*  AJVQ * AJHQ.  VQ.  CO.  ALSO 

COMMON  /ECS1/  AJVQ 
COMMON  /ECS2/  ajmO 
COMMON  /ECS3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECS5/  ALOO 


The  LOOPS  over  STATEMENT  8 CHANGE  THE  VERTICAL  CURRENT  AJV  FOR 
ALL  Tn£  MESH  POINTS  IN  THE  TRANSMISSION  PLATE. 


PLATE283 

PLATE28* 

PLATE285 

PLATE296 

PLATE287 

PLATE288 

PLATE289 

PLATE290 

PLATE291 

PLATE292 

PLATE293 

PLATE29* 

PLATE295 

PLATE296 

PLATE297 

PLATE293 

PLATE299 

PLATE300 

PLATE301 

PLATE302 

PLATE303 

PLATE30* 

PLATE305 


RESI2  aRESIS/2. 

OCLTIal ,/OELT 

00  2 Jal.JNCT 

CALL  REAOEC  (AJV(I).AJVQU.J).INCX) 

CALL  REAOEC  (VII). vQll.J), INCX) 

CALL  REAOEC  (AL(  I)  .ALQQll.J)  .INCX) 

00  1 I ■? ♦ I NCX 

AJVU)aIV»U-V«I-l>*AJViI)a»AL«I>*0ELTI-RESI2))/(ALII>*OELTI*RESI2PLATE306 
1*  - PLATE307 

CONTINUE  PLATE303 

AJV (1) a* AJV 12) PLATE309 

CALL  UPITEC  (AJV  ( 1 ) .AJVQ(I.J) . I NCX  > PLATE310 

CONTINUE  PLATE3II 

PLATE312 

The  LOOPS  OVER  STATEMENT  10  CHANGE  The  HORIZONTAL  CURRENT  AJH  FOR  PLATE313 
ALL  TnE  M£SH  POINTS  IN  Th Z TRANSMISSION  PLATE.  PLATE31* 

PLATE315 

00  * Ja2, INCT  PLATE316 

CALL  REAOEC  (AJH( 1 t .AJHO( I. J) . 1NCX)  PLATE317 

CALL  REAOEC  (V(II .VQ(|.J).INCXI  PLATE313 

CALL  REAOEC  IVld)*VQd«J-l)«INCX>  PLATE31 9 

CALL  REAOEC  I AL ( 1 ) < ALGU ( 1 « J) . I NCX)  PLATE320 

00  3 lal.INCX  PLATE321 

AJH( 1 1 a (V < I) »V1 ( II »AJh( I)  •(  AL( I) *UELT I “RES 1 2) ) / < AL ( l) *0ELTI*RESI2>  RLATE322 
CONTINUE  PLATE323 

CALL  KRITEC  lAJMll.AJHQ(l.J)tlNCX)  PLATE32* 

CONTINUE  PLATE325 


53 


t 


I 


I 


CALL  REAOEC  ( AJH( 1 ) 9 AJHQ (1*2) * INCA) 

C 

c The  loop  over  statement  ii  is  a reflected  boundary  condition 

C CAUSEO  BY  SYMMETRY. 

c 

00  S I»1 * INCX 

5 AJH(I>*-AJH(I> 

CALL  MRITEC  IAJHI 1) 9 AJHQ ( 1.1) • INCA) 

c 

C the  loops  over  statement  13  Change  The  VOLTAGE  v for  all  The  MESH 
C POINTS  IN  The  TRANSMISSION  PLATE. 

c 

DO  7 J«I.INCY 

CALL  READEC  IV ( 1) .VQ ( 1 . J> . INC*) 

CALL  REAOEC  UJHU>.AJrlO<l,J>.lNCAl> 

CALL  READEC  IAJV ( 1 ) . AJVO< 1. J> . INCXl) 

CALL  READEC  (C( 1> .CO M . J) * INCX) 

CALL  REAOEC  (VI  I )) .AJHU(1.J»1) .INCA1) 

00  6 Is I • INCX 

V 1 1 ) aV ( I ) *OELT*(AJV 1 1 ♦ 1 ) -A JV 1 1 ) • V l ( I ) - AJH ( I ) ) /C  < I ) 

6 CONTINUE 

CALL  VP I TEC  (V(I)»V0(1»J)* INCA) 

7 CONTINUE 

C WHEN  NSETUP  EQUALS  A or  when  NSETUP  EQUALS  5 

C Th«  LOOPS  OVER  STATEMENTS  1A  ANO  16.  17  AND  19  00  A REFLECTION 
C ACROSS  A AS  DEGREE  ANGLE  IN  Th€  CENTER  OF  ThE  CROSSED 

C TRANSMISSION  PLATE  70  IMPOSE  » FOLD  SYMMETRY. 

C 

IF  INSCTUP.NC.S.ANO.NSETUP.NE.A)  GO  TO  1A 
IF  » NSETUP. NE. 5. ANO. NSETUP. NE. A I GO  TO  1A 
00  • J»I.INCY1 

B CALL  REAOEC  <»B!l»J)»VQll»J>» INCY1) 

00  l«  LM»l * INCY 
J»INCY1-LM*1 
00  9 I at . J 

9 MSdf  JlaMBIJfl) 

CALL  aRITEC~Na8( 1 1 J) 9 Vd( It Jl » J) 

10  CONTINUE 

00  II  Jat.INCYl 

CALL  REAOEC  (aBIlfJ) 9 AJHO ( 1 9 J) 1 INCY 1 ) 

11  CALL  REAOEC  ( * A ( 1 9 J) 9 AJVO  < 1 * J) 9 INCY1 ) 

00  13  LMal.Is'CY 

J»INCY1-LM»1 
00  12  lal.J 
MS ( If J) aMA (Jf I ) 

12  MAdfJtaMBUf  1) 

CALL  MRITEC  (MB  ( 1 • J)  f AJrtQ  ( 1 9 J)  9 J) 

CALL  MRITEC  IMA  I If J) 9 AJVOi 1 9 J) 9 J) 

13  CONTINUE 

1A  CONTINUE 

IF  (NSETUP. NE. 2)  GO  TO  16 
CALL  REAOEC  (VdlfVOd.D.INCX) 

CALL  REAOEC  (VKII  fVOdfiNCTDtINCXt 


54 


; — ~ ■ 

'■.-•■WAK . 


PLATE326 
PLATE327 
PLATE323 
PLATE329 
PLATE330 
PLATE331 
PLATE332 
PLATE333 
PLATE33A 
PLATE33S 
PLATE336 
PLATE337 
PLATE338 
PLATE339 
PLATE3A0 
PLATE3AI 
PLATE3A2 
PLATE3A3 
PLATE3AA 
PLATE3A5 
PLATE3A6 
PLATE3A7 
PLATE3A6 
PLATE3A9 
PLATE350 
PLATE351 
PLATE352 
PLATE3S3 
PLATE35A 
PLATE35S 
PLATE356 
PLATE357 
PLATE356 
PLATE359 
PLATE36J 
PLATE361 
PLATE362 
PLATE363 
PLATE36A 
PLATE365 
PLATE366 
PLATE367 
PLATE368 
PLATE369 
PLATE370 
PLATE371 
PLATE372 
PLATE373 
PLATE37A 
PLATE3TS 
PLATE 376 
PLATE377 
PLATE37B 
PLATE379 


V 


I 


i 


15 

16 


00  15  l«ltINCX 
VU)«Vl(!>»<tf<I)»Vl<H  1/2.0 


CONTINUE 
CALL  naiTEC 
CALL  «BITEC 
CONTINUE 


(V(l)tVQ(ltll* INCX) 
<V1U> t VO  1 1 * I NCT 1 ) * I NCX ) 


return 


ENO 


PLATE380 

PLATE381 

PLATE382 

PLATE383 

PLATE38* 

PLATE385 

PLATE386 

PLATE38T 


c 

c 

c 

c 

c 

c 


SUBROUTINE  MISCELL 


COMMON  X. 

Rl. 

VO. 

INCA. 

MOIST  1 • 

OISTO. 

RESIS. 

AJV (26) 

1 

T. 

R2« 

CT. 

INCT. 

N01ST2* 

IMCXO. 

KSEMS*. 

AJH (26) 

K« 

R3. 

ALT. 

AINCX. 

N0IST3. 

TSTOP. 

nsetup. 

V (26) 

T. 

R*. 

NT. 

AINCT. 

MOIST*. 

TEDIT. 

NSETUP2. 

C (26) 

AK. 

RAO. 

ALN. 

INCA  1 . 

OISTI. 

UELEO. 

MSETUP3. 

VI (26) 

OCLT. 

AL1. 

ALO. 

INCY1. 

0IST2. 

TE0IT2. 

MSETUP*. 

AL (26) 

KSCN. 

AL2. 

ANU. 

OELX. 

0IST3. 

MP0SA1. 

NSETUP5. 

•A(26.26) . 

NCAP. 

N100. 

NZAP  . 

OELT. 

OIST*. 

MP0SA2. 

moot . 

#8(26.26) . 

At ( 180) . Yl(100>. 

AT ( 1)  • 

CTRAC(I) 

• CTR 1(1) 

. CTR2U) 

01 MENS  ION  AJVQ (26*26) • AJmO ( 26*26) • VQ<26.26).  CQ(26*26>* 

1 ALQO (26*26) 

LEVEL  3*  AJVQ*  AJriO.  VO*  CO*  ALQQ 

COMMON  /ECSl/  AJVQ 
COMMON  /ECS2/  AJMO 
COMMON  /ECS3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECSS/  ALQO 

this  subroutine  sets  up  the  rlate  problem. 

REAO  STATEMENTS  ANO  SNUFFLING  of  REAL  to  integer  variables. 


PLATE388 

PLATE389 

PLATE390 

PLATE391 

PLATE392 

PLATE393 

PLATE39* 

PLATE395 

PLATE396 

PLATE397 

PLATE393 

PLATE399 

PLATE*00 

plateaoi 

PLATEA02 

PLATEA03 

PLATEaO* 

PLATEAOS 

PLATEA06 

PLATEAOT 

PLATE*06 

PLATE*09 

plateaio 
PLATEAU 
PLATE*! 2 
PLATE*13 
PLATEAJA 
PLATEA15 
PLATEA16 


PRINT  S 

DCLZ*6.«015 

ANTalOOOS. 

INCA  *26-1 
INCT  *26  - 1 

REAO  7*  PI «R2*RA0*R3*ER * SETUP *SETUP2» SETUP 3 
REAO  T.  X.T«VO*CT. ALT. CLENGTH. FACTOR. FACT0R2 
REAO  T.  TSTOP.OELZZ.OISTO.OISU.OIST2.01ST3.0ISTA 
IF  (0CL2Z.NC.0.0)  GCLZ*OCLZZ 

AINCA»INCA 

AINCY»INCT 

AZ-AINCY/T 

NSETUP-SETUP 

NSETUP2*SETUP2 

NSETUPJaSETUP3 

IF  (NSCTUP.EO.D  PRINT  22 

IF  INSETUP.EQ.2)  PRINT  23 

IF  (NSCTUP.CQ.3)  P«tNT  2* 

IF  (NSCTUP.Ca.*)  PRINT  2S 
IF  (NSCTUP.CQ.St  PRINT  2* 

PRINT  21.  Rl .«2. RAO. R3.CR. SETUP. SETUP2. SETUPS 
PRINT  21.  A.Y.VO.CT.ALT.CLENOTM.FACTOR.FACTOR2 
PRINT  21.  TST0P.0CL2Z.0IST8.0ISU.01ST2.0IST3.01ST* 

RCS1Sa|.67C-8* 

TCOIT*0.6*C-86 


PLATEai 7 
PLATEAU 
PLATEA19 
PLATEA20 
PLATEA21 
PLATEA22 
PLATEA23 
PLATEA2* 
PLATE A25 
PLATEA26 
PLATEA27 
PLATEA28 
PLATEA29 
plate *3o 
PLATEA31 
PLATEA32 
PLATEA33 
PLATEA3* 
PLATEA3S 
PLATEA36 
PLATEA37 
PLATEA38 
PLATEA39 
PLATE** 3 
PLATE ** 1 
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I 

» 

I 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


DCLEDaO. 021-06 

TC0IT2-0. 011-06 

RCAO  !•  A. 0.0 

IF  U.NE.0.0)  TEOITaA 

IF  (B.NE.O.O)  OELEJaS 

IF  (O.NE.O.O)  TEDIT2-0 

PRINT  2t  A. 0.0 

PRINT  3.  TEOIT. OILED. TE0IT2 

PRINT  14.  R 1 .R2.RAO.R3 

PRINT  15.  NSETUP. NSETUP2.NSETUP3 

PRINT  16t  X.V. FACTOR. FACT0R2.ER 

PRINT  17.  VO.CT .ALT .CLENGTH 

PRINT  18.  TSTOP.ANT.AInCX.AINCY 

PRINT  19.  0IST1.0IST2.0IST3.0IST-.0IST0 

NSETUP  a 1 P900UCES  A SQUARE  TEST  PROBLEM. 

NSETUP  a 2 PROOUCES  A PARALLEL  PLATE  TRANSMISSION  LINE  PROBLEM . 
NSETUP  a 3 PROOUCES  A RECTANGULAR  PROBLEM  *IT*  CAPACITOR  BANKS 
ON  TaO  SIDES. 

NSETUP  a 4.  PROOUCES  A RECTANGULAR  PROBLEM  «ITM  CAPACITOR  BANKS 
ON  ALL  FOUR  SIDES. 

NSETUP  a 5 PROOUCES  Th£  CROSSEO  PLATE  GEOMETRY. 


INCX  ANO  INCt  ARE  ThE  NUMBER  OF  X ANO  Y INCREMENTS  RESPECTIVELY. 

INCXIaINCX*! 

INCY laINCY ♦ 1 
INCXOalNCX-1 

NT  IS  The  number  OF  TIME  INCREMENTS. 

NTaANT 

PI*3. 1-15926S35 

AMUO  IS  The  PERMEABILITY  OF  FREE  SPACE. 

EPSI  IS  T*E  PERMITTIVITY  OFMYLAR. 

AMUOa4.0*PI*4.0E-07 

EPSI0«1.0/ ( AMUO* I2.99T9E-08) **2> 

EPSIa£R*EPSIO 

PRINT  2«o  AMUO. EPSIO. EPSI 

CEL*  ANO  OELY  ARE  THE  DIMENSIONS  OF  ONE  CELL.  CEL*  ANO  OELY 

smoulo  be  approximately  eoual.  del 2 is  the  plate  separation. 

CELXaX/ < AINCX-1 .0) 

OELYaY/ I AINCY— 1 .01 

NPOSXlal <0ISTl«0IST2)/2.0t/0ELX 

NP0SX2«<  <OIST3«OIST4)/2.0)/OEl* 

PRINT  4.  NPOSX 1 .NP0SA2 

NOISTI  ANO  N01ST2  ARE  THE  X COORDINATE  BOUNDARIES  OF  ONE 
CAPACITOR  BANK  MOOULE  ON  TmE  SIOE  OF  THE  CROSSED  TRANSMISSION 


PLATE442 
PLAIE443 
PLATE444 
PLATE445 
PLATE446 
PLATE447 
PLATE44B 
PLATE449 
PLATE-53 
PLATE451 
PLATE452 
PLATE453 
PLATE454 
P LATE455 
PLATE456 
PLATE457 
PLATE458 
PLATE459 
PLATE-60 
PLATE-61 
PLATE-62 
PLATE463 
PLATE464 
PLATE-65 
PLATE-66 
PLATE-67 
PLATE468 
PLATE-69 
PLATE-70 
PLATE471 
PLATE- 72 
PLATE473 
PLATE-74 
PLATE-75 
PLATE-76 
PLATE-77 
PLATE-78 
PLATE-79 
PLATE-83 
PLATE-81 
PLATE-82 
PLATE-83 
PLATE-8- 
PLATE-85 
PLATE-e6 
PLATE-87 
PLATE-88 
PLATE-89 
PLATE-90 
PLATE-91 
PLATE-92 
PLATE-93 
PLATE-9- 
PLATE-95 
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i 


! 

0 

c 

PLATe. 

PLATEA9* 

- 

c 

N01ST3  ANO  N0IST4  ARC  THE  X COORDINATE  BOUNDARIES  OF  THE  SECONO 

PLATEA97 

c 

MODULE  ON  The  SIDE  OF  TmE  TRANSMISSION  PLATE. 

PLATEA9B 

c 

PLATE499 

NOISTUIOIST1/OCLO*! 

PLATE500 

N0tST2>(0IST2/0ELX)»l 

PLATE501 

NO 1 ST3» ( D t ST  3/ DEL A ) •! 

PLATES02 

NOISTAafOISTA/OELA) ♦ ! 

PLATE503 

c 

PLATE50A 

c 

AL2  ANO  C2  ARE  ThE  1N0UCTANCE  ANO  CAPACITANCE  IN  THE  NORMAL  PART 

PLATE505 

c 

OF  A TRANSMISSION  PLATE. 

PLATE506 

c 

PLATES07 

AL2*AMU0*0ELY*0ELZ/DELX 

PLATE508 

C2»£PS I ‘OELA  «OEL Y/OELZ 

PLATE509 

c 

PLATE510 

l 1 e 

AM  ANO  Cl  ARE  THE  INOUCTANCE  ANO  CAPACITANCE  OF  EACH  CAPACITOR 

PLATE51 1 

c 

ANO  SWITCH  ASSEMBLY  ASSUMING  ThE  NUMBER  OF  CAPACITORS  EOUALS  ThE 

PLATE512 

c 

NUMBER  OF  * MESH  POINTS  WITHIN  THE  BOUNDARIES  OF  A CAPACITOR  BANK 

PLATE513 

c 

MODULE. 

PLATE51 A 

c 

PLATE515 

NC APmCL  ENG  T H/ DEL Y 

PLATE516 

IF  (NCAP.E9.0)  NCAP»INCY 

PLATE517 

IF  (NCAP.GT.INCYl  ncap»incy 

PLATE519 

CI*CT/(NCAP-l.l 

PLATE519 

ali«altxncap-i.o> 

PLATE520 

c 

PLATE521 

c 

ALA  ANO  C4  ARE  The  INDUCTANCE  AND  CAPACITANCE  IN  THE  LOAD  WHERE 

PLATE522 

c 

C4  IS  SET  ARBITRARILY  LARGE. 

PLATES23 

c 

PLATE52A 

* 1 

C4«c2*i«ooooe. 

PLATE525 

PLATE52S 

AL4-AL2 

PLATE527 

c 

PLATES28 

c 

CELT  IS  Th£  TIME  STEP 

PLATE529 

c 

PLATE530 

DELT »P I / 20 . »SOR T « C2» AL2 ) 

PLATC531 

OELT»0£LT»2.0 

PLATES32 

PRINT  B 

PLATE533 

PRINT  B.  AlNCA » AINCY * ANT • X » Y 

PLATE53* 

PRINT  9.  VO.CT.AtT 

PLATES3S 

c 

PLATE538 

c 

A ANO  T APE  THE  DIMENSIONS  OF  THE  TRANSMISSION  PLATE  (METERS). 

PLATE537 

c 

oelz  is  the  thickness  of  the  mylar. 

PLATE538 

c 

PLATES39 

PRINT  Ilf  C2f AL2 

PLATE5A0 

PRINT  12t  C1.AL1 

PLATESAI 

PRINT  13f  C4.AL4 

PLATEShZ 

PRINT  I0f  OELT tOCLAf DELYtOELZt RI 

PLATE543 

PRINT  6 

PLATESAA 

CALL  MCSM  (Cl tC2fC3fCA#AL3f ALA. FACTOR. FACTOR 1 fFACTOR2 I 

PLATE5A5 

CALL  L0A0ISETUP3) 

PLATESAA 

CALL  WfDOC  (Ci.C2.C3.CA.AL3. ALA. FACTOR. FACT0R1 .FACT0R2) 

PLATE5A7 

V-ve  <,"• 

CALL  FRAME 

PLATESAA 

CALL  PRINT  1 (1) 

PLATES49 

• 

53 

% 

[ ' 

• -'..'U-mi  «-■  . m . 

✓ 

30 

31 

SO 

32 


C 

c 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 
1* 

15 

16 

17 

18 
IS 

20 

21 

22 

23 


CALL  PR INTI  (2) 

CALL  PRINT  1 (3) 

PRINT  6 

FORMAT<1#X,*OIAOONAL  CAPACITIVE  ELEMENTS  ARE*./) 
PRINT  30 
FORMAT  1 10E 10.2) 

INCTO»INCr-l 

00  50  1*1* INCYQ 

CALL  REA0EC1X1 ( I ) »CQ (I i I ) « 1 ) 

CALL  PEAOEC (T1 ( I ) « ALQO ( It  I) 1 1) 

CONTINUE 

PRINT  31 . <X1(I) .1-1.1NCY0) 

FORMAT <10X. *01  AGONAL  INDUCTIVE  ELEMENTS  ARE**/) 
PRINT  32 

PRINT  31. <T1CI).I«1.INCT0) 

RETURN 


PLATE550 
PLATE5S1 
PLATE552 
PLATE5S3 
PLATE554 
PLATE5SS 
PLATE556 
PLATE557 
PLATE5S3 
PLATE559 
PLATE560 
PLATE561 
PLATE562 
PLATE563 
PLATE56- 
PLATE565 
PLATE566 
PLATE567 
PLATE568 
PLATE569 
TE0IT2  «PLATES70 
PLATE571 
PLATES72 
PLATES73 
PLATE574 
PLATE57S 
> .E14.6.PLATE576 
PLATE577 
.1PEIA.PLATES78 

. _ PLATE579 

FORMAT  (3A.7M0ELT  • .IPC14.6.8H  OtLX  * .E14.6.3H  OELY  ■ .E14.6.8H  PLATE580 
10ELZ  a .E14.6.24H  FRACTIONAL  HOLE  SIZE  ■ ,£14. 6./)  PLATES81 

FORMAT  I3X.SSHCMARACTERISTICS  OF  Th£  TRANSMISSION  PLATE  ELEMENTS  APLATE582 
IRE  . 1PE1A.6.9H  FARACS  ,E|4.6,9H  HENRIES../)  PLATE583 

FORMAT  (AX.72HCMARACTEPISTICS  OF  EACH  mESH  POINT  THAT  SERVES  AS  A PLATE53* 
lBREAKOONN-  POINT  arc  . 1PC14.6.9H  FARADS  .El4.6.9n  HENRIES../)  PLATE585 

FORMAT  (SX.61HCHARACTERISTICS  OF  EACH  MESH  POINT  THAT  SERVES  AS  A PLATES86 
1SINK  are  »IP€1«».6.9M  FARAOS  .E14.6.9H  HENRIES../)  PLATES87 

FORMAT  <10X.ShR1  a .1PE1A.6.6H  R2  a ,1PE14.6.7h  RAO  a .1PE14.6.6H  PLATE586 
1R3  « . IPElArO./)  PLATES89 

FORMAT  ( 19X.9m.NSE TUP  a ,I3,UH  NSETUP2  a ,13, Hh  NSETUP3  a .13,/)  PLATE590 
FORMAT  <1»X,4HX  a , 1PE1A.6.5H  Y a ,1PE14.6,10H  FACTOR  a , 1PE14.6. 1PLATE591 
1IM  FACT0P2  ■ .1PEN.6.6H  ER  a .1PE14.6./1  PLATE5S2 

FORMAT  •SX.SHVO  a , 1PE1A.6,6h  CT  a ,1?£14.6.7h  ALT  a , 1PE14.6. 10HCPLATE593 
1LEN6TH  a , 1PC14.6./)  PLATE594 

FORMAT  (SX,8hTST0P  a , IPE 14.6, 7H  ANT  a , 1P£14.6*9h  AINCX  a , 1PE14.PLATE595 
I6.9H  AINCY  a , 1P£1*,6./) 

FORMAT  (SX«8hOIST 1 a ,lPtl4.6.9H  01ST2  a , 1PE14.6.9H  0IST3 
14.6.9m  0IST4  a ,1P£1».6.9h  OISTO  ■ .l°£l4.6./> 

FORMAT  (9X.7HAMU0  a ,lt»El4.6,9H  EPSIO  a ,1P£14.6.3H  EPSI 


FORMAT  (8F10.2) 

FORMAT  (13X.2P6E12.1,/,/) 

FORMAT  ( 19X.6HTE0IT  a .2PE12.1.10H  DELED  a .2PE12.1.11H 
1 .2PE12. 1 ././) 

FORMAT  (17X.9HNP0SX1  a ,15,19m  NP0SA2  a ,IS,/) 

FORMAT  <1h1> 

FORMAT  (1H0) 

FORMAT  (8F19.2) 

FORMAT  ( IX ,7hINCX  a .1PE14.6.8H  INCY  a .£14,6. 10H  4 
12SH  X, Y PLATE  DIMENSIONS  a .£14.6. 3H  X .£14.6./) 

FORMAT  (2X.45nCHARACTERISTICS  OF  ONE-HALF  OF  A MOOULE  ARE 
16. 3H  VOLTS  .£14.6. 13H  FARAOS  ANO  .£14.6. 8h  HENRIES./) 


I NCR 


16./1 

FORMAT 


(I4X,1P8E12.2./) 


FORMAT  <DX,37mSSSSS 
FORMAT  (10XO8HSSS66 


SYMMETRY  TEST  PROBLEM 
P.P.  TRANSMISSION  line 


SSSSS./) 

tmt,/) 


PLATES96 
.1PE1PLATE597 
PLATE598 
.1PE14.PLATE599 
PLATE600 
PLATE601 
PLATE602 
PLATE603 
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f 


format 

format 

format 

ENO 


<14Xt«*HSSSSS  T«0  SIDED  TRANSMISSION  PLATE  SSSSS*/) 

1 1 9X»*5HtSSiS  FOUR  SIOEO  TRANSMISSION  PLATE  SSSSS'/) 
mXOSMSSSSS  CROSSED  PLATE*  l MJ  SSSSS'/) 


PLATE60A 

PLATE605 

PLATE606 

PLATE607 


SUBROUTINE  MESH<C1 .C2.C3.C4.AL3.AL4.FACT0R.FACT0R1 .FACTOR2) 


PLATE608 

PLATE609 


COMMON  X. 

»1. 

vo. 

INCX. 

NOISTl. 

OISTO. 

RESIS. 

AJV (26) . 

PLATE610 

1 

T. 

R2. 

CT. 

INCT. 

N0IST2. 

INCXt. 

KSENS*. 

AJH(2b>  . 

PLATE61 1 

2 

K. 

R3. 

ALT. 

A INCX. 

N0IST3. 

TSTOP. 

NSETUP. 

V (26) • 

PLATE612 

3 

T. 

R4. 

NT* 

AINCT. 

N0IST4* 

TEOIT. 

NSETUP2. 

C (2b)  • 

PLATE613 

4 

AK. 

RAO. 

ALN. 

INCX I . 

0IST1. 

oeleo. 

NSETUP3. 

VI (26). 

PLATE614 

s 

OELT. 

AL). 

ALQ. 

INCTI. 

UIST2. 

TE0IT2. 

NSETUP4, 

AL (26) . 

PLATE61S 

6 

KSEN. 

AL2. 

ANU. 

OELX. 

0IST3. 

NPOSX 1 • 

NSETUPS. 

■A (26.26) . 

PLATE616 

7 

NCAP. 

N100. 

NZAP. 

CELT. 

0IST4. 

NP0SX2. 

1 1 1 ( 10>  • 

*8(26.26)  . 

PLATE617 

8 

XI (100) 

• ri(ioo). 

AT ( 1) . ' 

CTRAC  ( 1 ) 

. CTRI(l) 

. CTR2( 1 ) 

PLATE6I8 

OIMENSION  AJVQ(26.26>  » AJM0(26.26> » VO(26.26).  CU126.26). 
I ALOO (26.26) 

LEVEL  3.  AJVQ*  AjnU.  vg.  CO.  ALQO 

COMMON  /EC SI/  AJVQ 
COMMON  /ECS2/  AJMU 
COMMON  /ECS3/  va 

common  /£cs*/  ca 

COMMON  /ECSS/  ALOO 


this  loop  initiallizes  capacitance  and  voltage  throughout  the 

MESH.  , 

ZERO *6.0 

00  I Jal.INCTl 

oo  i i-i.iNCxi 

CALL  WRITEC  (C2.C0( I t J> . 1) 

CALL  HRITEC  (ZERO. VO ( I. J) • 1 ) 

CALL  ARITEC  (ZERO. AJV0( I • J) • 1 ) 

CALL  WRITEC  (ZERO* AJho< I • J) • 1) 

CALL  WRItEC  ( AL2. ALOO( l. J) • 1) 

continue 

IF  (NSETUR.EQ.2)  GO  TO  7 ~ 

This  SECTION  SETS  UP  A TEST  PROBLEM  TO  OETERMInE  AZIMUTHAL 
ASSTMMETRV  when  NSETUP  ■ 1.  THIS  IS  THE  ASSTMmETHY  WHICH  IS 
CAUSEO  BT  The  SOU ARE  h£Sh. 

C2FAC*C2/F ACTOR 

ALFAC*AL2«FACT0R 

IF  (NSETUP.NE.l)  00  TO  3 

00  2 J»1.1NCY1 

00  2 1*1. INCA  1 

AX*l*OELX 

AT*J«OELT 

R AO -SORT ( AX*AX. AT* AT ) 

IF  (RA0.0T.R3)  CAU.  •RUCC  (Cl .CQ ( I . J) . 1) 

IF  IRA0.GT.R3)  CALL  RR1TCC  (VO.VO(I.J).l) 

ALMMM*ALl/S. 

CT10*CTM0. 


PLATE619 

PLATE620 

PLATE621 

PLATE622 

PLATE623 

PLATE62* 

PLATE62S 

PLATE626 

PLATE627 

PLATE626 

PLATE629 

PLATE630 

PLATE631 

PLATE632 

PLATE633 

PLATE63* 

PLATE635 

PLATE636 

PLATE637 

PLATE638 

PLATE639 

PLATE640 

PLATE641 

PLATE642 

PLATE643 

PLATE644 

PLATE64S 
PLATE646 
PL  A T E 64  7 
PLATE646 
PLATE649 
PLATE650 
PLATE651 
PLATE6S2 
PLATE653 
PLATE6S*. 
PLATE6S5 
PLATE6S6 
PLATE657 
PLATE658 
PLATE659 
PLATE660 
PLATE661 


uuuuuu  <f  uuuu  in  ouuu 


F(RaD.GT.R3-3.*0ELY>  CALL  MR 1 TEC < ALMMM* ALQO (l.J)tl) 
F (RA0.LT.R2)  CALL  MRITEC  ( ALFAC* ALOQ ( l • J) • 1) 


(RA0.LT.R2) 

(RA0.LT.R1) 


CALL 

CALL 

CALL 


MRITEC 
MR  I TEC 
MR1TEC 


(C2FAC.CUU.J)  .1) 
(cnotcodtj).i) 
(AL2.AL00U.J) . 1) 


r (RA0.LT.R1) 

2 CONTINUE 
GO  TO  10 

3 CONTINUE 
LN»INCT/2*1 

the  loops  over  inoex  a set  up  a region  of  greater  plate 

JUST  OUTSIOE  The  hOLE  IN  TmC  CENTER  MmICH  REPRESENTS  THE  LOAO. 
The  variable  factor  IS  the  ratio  of  THE  nEM  SEPARATION  TO  the 

0L0  SEPARATION. 

00  A I«l.LN 
00  A J»1.LN 
AX«I*0CLX 
AT»J*OELT 

RAaSORT(AX*AX»AY*AT) 

CT10«CT*10. 

IF  (RA.LT.R2)  CALL  MRITEC 
IF  (RA.LT.R2)  CALL  MRITEC 
IF  (RA.LT.Rl)  CALL  MRITEC 
IF  (RA.LT.Rl)  CALL  MRITEC 
CONTINUE 
I-INCX1 


PLATE662 
PLATE663 

PLATES6A 
PLATE665 
PLATE666 
PLATE66T 
PLATE666 
PLATE669 
PLATE670 
PLATE671 
SEPARATI0NPLATE672 


(C2FAC.C0( I. J) . 1) 

( ALFAC • ALOQ ( I > J) » 1 ) 
(CA.CU(I.J) .1) 

( ALFAC. ALQOd. J)  .1) 


The  LOOP  over  5 CHARGES 
TRANSMISSION  PLATE. 


The  CAPACITORS  AT  the  ENO  OF  The 


00  S Jal.NCAP 

CALL  MRITEC  (Cl*CO(I.J).ll 

CALL  MRITEC  (V0«V0(I.J>'1) 

CALL  MRITEC  (ALl*ALQO(l.J> .1) 

CONTINUE  — 

IF  (NSETUP.NE.S)  GO  TO  10 

The  LOOP  OVER  7 CHARGES  The  TmO  CAPACITOR  BANK  MODULES  ON  THE 
SIOC  OF  THE  TRANSMISSION  PLATE. 


JMINCYI 

00  6 lal.INCX 

001ST»I»0€LX 

IF  ( 00 1 ST .GE. 01 ST  1 .ANO.QOIST .LE.0IST2) 
(00IST.GE.0IST1 .ANO.GOIST .LE.0IST2) 
(00 1ST .GE.OIST 1, ANO.OOIST .LE.01ST2) 


IF 
IF 
11) 

IF 
IF 
IF 
11) 
CONTINUE 
GO  TO  10 


(OOIST.GE.OIST3.ANO.QOIST.LE.OIST*) 

(00IST.GE.0IST3.AN0.001ST.LE.0IST*) 

(001ST.GE.DIST3.AN0.g0IST.LE.0lST*) 


CALL  MRITEC 
CALL  MRITEC 
CALL  MRITEC 

CALL  MRITEC 
CALL  MRITEC 
CALL  MRITEC 


PLATE673 
PLATE67*. 
PLATES75 
PLATE67S 
PLATE677 
PLATES78 
PLATE679 
PLATE680 
PLATE681 
PLATE682 
PLATE683 
PLATE68** 
PLATE68S 
PLATE686 
PLATE687 
PLATE688 
PLATE689 
PLATE69Q 
PLATE691 
PLATE692 
PLATE693 
PLATE69* 
PLATE695 
PLATE696 
PLATE697 
PLATE698 
PLATE699 
PLATE700 
PLATE701 
PLATE702 
PLATE703 
PLATE70* 
PLATE705 
(Cl.CO(l«J) *1)  PLATE706 
(V0*V0(I*J>«1>  PLATE707 
(ALliALOO(ItJ) .PLATE708 
PLATE709 
(CltCQ(I.J) . 1)  PLATE710 
(VO.VQ ( I . J) • 1 ) PLATE711 
( AL 1 • ALOC ( I • J ) • PL A T E 7 1 2 
PLATE713 
PLATE71* 
PLATE715 


62 


•jfmon 


WHEN  NSETUP-2  THIS  SETS  UP  A PARALLEL  PLATE  TRANSMISSION  LINE 
PROBLEM  FOR  COMPARISON  TO  AN  ANALYTICAL  SOLUTION. 

CONTINUE 
1*1 

00  B J»l  • INCTI 
CALL  WRITEC 
• CONTINUE 

I-INCX1 

INQ-OISTI/OELTd 

INQOaOIST2/OCLT»l 

nu*«inoo-inqm 

CtaCT/NUM 

AL1«ALT»NUM 

00  9 J*I NO « INOQ 

CALL  MRITEC  <Cl«CJ<I*J)tl> 

CALL  WRITEC  (VO.YO(I.J) .1) 

9 CONTINUE 

10  CONTINUE 
return 
ENO 


PLATE7I6 

PLATE717 

PLATE718 

PLATE7I9 

PLATE720 

PLATE721 

PLATE722 

PLATE723 

PLATE72A 

PLATE725 

PLATE726 

PLATE727 

PLATE729 

PLATE729 

PLATE730 

PLATE731 

PLATE732 

PLATE733 

PLATE73N 

PLATE73S 

PLATE736 

PLATE737 
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nn  nnn 


I 


SUBROUTINE  LOAD  (SETUP3) 


PLATE738 

PLATE739 


COMMON  X. 

Rl. 

VO. 

INCX. 

NO I ST l . 

OISTO. 

RESIS, 

AJV (26) . 

PLATE7A0 

1 

Y. 

R2. 

CT. 

INCY. 

NDIST2, 

INCXO. 

KSENSA. 

AJHI26) . 

PLATE7A1 

2 

K . 

R3. 

ALT  . 

AINCX. 

N0IST3. 

TSTOP. 

nsetup. 

V (26) . 

PLATE7A2 

3 

T. 

PA. 

NT. 

AINCY. 

NOISTa. 

TEDIT. 

NSETUP2, 

C (26)  . 

PLATE7A3 

* 

AK. 

RAO. 

ALN. 

INCX 1 . 

0IST  | . 

DELED. 

NSETUP3. 

VI (26) . 

PLATE7A*. 

5 

OELT « 

AL1. 

ALQ. 

INCY  I . 

0IST2. 

TEDIT2. 

NSETUPA, 

AL (26) . 

PLATE7A5 

6 

KSEN. 

AL2. 

ANU. 

OELX. 

0IST3. 

NP0SX1. 

NSETUP5, 

•A ( 26. 26) • 

PLATE7A6 

7 

NCAP  « 

N100. 

NZAP. 

OELY. 

OISTA. 

NP0SX2. 

1 1 1 < 10)  • 

•8(26.26) . 

plate  7a7 

8 

Xl(lOO) 

. YK100). 

AT  ( 1 ) . 

CTRAC(l) 

. CTRUl) 

. CTR2U) 

PLATE7A8 

DIMENSION  AJVQ(26.26).  AJHQ<26.26>,  VQI26.26).  C3«26.26>. 
1 ALQQ (26.26) 

LEVEL  3.  AJVQ.  AJHQ.  VQ.  CQ.  ALQQ 

COMMON  /ECSl/  AJVQ 
COMMON  /ECS2/  AJHQ 
COMMON  /ECS 3/  VQ 
COMMON  /ECS*/  CQ 
COMMON  /ECS5/  ALQQ 


TmIS  SECTION  INCREASES  The  INOUCTANCE  IN  The  LOAO  TO  SIMULATE  A 
SHIVA  LOAO  INSTEAO  or  A SHORT  CIRCUIT.  (AhEN  SETUP3.NE.0.O. 
ALM.J)*SETUP3.) 

IP  ( SETUP 3. £3 .0.0)  <30  TO  2 
NNNa AMINl ( AlNCX • AINCT « 55. ) 

00  1 I ■ 1 . NNN 
00  1 Jat.NNN 
AX»I»OELX 
AY«J*OELT 

RAaSQRT ( AX»AX»AY»AY) 

ir  OA.LT.Rl*3.0»0ELX)  CALL  VRITEC  ( SETUP 3.  ALQQ  ( I . J)  . 1) 

CONTINUE 

CONTINUE 

RETURN 

ENO 


PLATE7A9 
PLATE750 
PLATE751 
PLATE7S2 
PLATE753 
PLATE75A 
PL  ATE  755 
PLATE756 
PLATE757 
PLATE7S8 
PLATE759 
PLATE760 
PLATE761 
PLATE762 
PLATE763 
PLATE76-. 
PLATE765 
PLATE766 
PLATE767 
PLATE768 
PLATE769 
PLATE770 
PLATE771 
PLATE772 
PLATE773 
PLATE77* 
PLATE775 
PLATE776 
PLATE777 


J 


c 

c 

c 

c 

c 


SUBROUTINE 

WEDGE 

(Cl* 

C2*C3»C4«al3*al4* factor* 

F ACT0R1 .FACT0R2) 

PLATE773 

PLATE779 

COMMON  X. 

Rl* 

VO* 

INCX. 

MOIST  1. 

OISTO. 

RESIS. 

AJV (26) * 

PLATE780 

I 

Y* 

R2. 

CT. 

INCY* 

N01ST2. 

1NCX0. 

KSENSW. 

AJH (26) • 

PLATE781 

2 

K * 

R3. 

ALT* 

AINCX. 

NO  1ST  3* 

TSTOP. 

NSETUP. 

V (26) » 

PLATE782 

3 

T* 

R4* 

NT* 

aincy. 

MOIST  4* 

TEOIT. 

NSETUP2. 

C (26) • 

PLATE783 

4 

AK* 

RAO* 

ALN* 

INCX1* 

01ST 1 » 

DELED* 

NSETUP3. 

VI (26) « 

PLATE78*. 

s 

DELT, 

AL1* 

ALO* 

INCY 1 • 

DIST2, 

TEDIT2, 

NSETUP* , 

AL (26) * 

PLATE78S 

6 

KSEN. 

AL2* 

ANU* 

DELX* 

0IST3. 

NP0SX1. 

NSETUP5. 

WA126.26) • 

PLATE786 

7 

NCAP  * 

NlOO* 

NZAP. 

OELY. 

0IST4. 

NP0SX2* 

I I l ( 10) • 

w6(26*26) • 

PLATE787 

8 

XI  MOO)  • Yl  ( 106)  . 

ATM). 

CTRAC  M) 

. CTR1M) 

• CTR2M) 

PLATE78R 

DIMENSION  AJVO (26*26)*  AJHQ<26*26>.  vg<26.26>.  CJ(26*26)* 

1 ALQO (26*26) 

LEVEL  3*  AJVQ,  AJhOt  VO*  CO*  ALQO 

COMMON  /ECS1/  AJVQ 
COMMON  /ECS2/  AjrtQ 
COMMON  /ECS3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECS5/  ALSO 

RK«I 

CAUL  SETUPX  (1> 

CONTINUE 

READ  6*  M*N*K*L*NQ 
IF  (K.EQ.O)  60  TO  3 
IF(NO.GT.O)  FACT0R2-NO 

•MEN  THIS  SECTION  IS  USEO  CARE  IS  NECESSARY  TO  INSURE  THAT  KK 

NEVER  EXCEEDS  THE  DIMENSIONED  SIZE  OF  XI  ANO  Y1  ARRAYS* 

if  that  ooes  happen*  some  of  the  common  block  will  be  wiped  out. 

PRINT  5*  M.N«K*L*NQ 
IF  (L.EQ.O)  L-INCX 
C2F2«C2/F ACT0R2 
AL2F2»AU2»FACT0R2 
00  2 I >K  »L  — 

00  2 J»H,N 

CALL  wRITEC  <C2F2*C0 ( I » J) • 1> 

CALL  WRITEC  ( AL2F2* ALOU ( I • J) » 1) 

CONTINUE 

XI(KKJ»r«0ELX 

tkkk»«m*oely 

xkkk*i>»k*OEH 

YICKK*1»»N*0ELY 
X1(KK*2>«L*0ELX 
Y1  (KK«2) «N#OELY 
XI (KKO)  «L#OELX 
Yl (KK*3I »M*OELX 
XI (KK*4I PX 1 IKK) 

Y1(KK»*)«Y1IKK) 

K«5 


P LATE7B9 
PLATE790 
PLATE791 
PLATE792 
PLATE793 
PLATE 79* 
PLATE795 
PLATE 706 
PLATE797 
PLATE793 
PLATE799 
PLATE800 
PLATE30I 
PLATE802 
PLATE303 
PLA7E60** 
PLATE80S 
PLATE806 
PLATE807 
PLATE808 
PLATES09 

plates i o 
PLATE8U 
PLATE812 
PLATE813 
PLATE31* 
PLATE315 
PLATE816 
PLATE317 
PLATE818 
PLATE3I9 
PLATE320 
PLATE821 
PLATE822 
PLATE323 
PLATE82* 
PLATE825 
PLATE326 
PLATE327 
PL ATE 826 
PLATE829 
PLATE830 
PLATE831 
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p 


call  setup* (2> 

GO  to  i 

3 CONTINUE 

return 


format  ( 8F 10.2) 
format  (2X*5I5) 
FORMAT  1612) 

ENO 


PLATC632 

PLATE833 

PLATES}* 

PLATES3S 

PLATE836 

PLATE637 

PLATES38 

PLATES39 

PLATE840 

PLATES*! 


SUBROUTING  SETUPX  <KLH> 


PLATE842 

PLATE843 


COMMON  X* 

«!• 

VO* 

INCX. 

N0I5T I * 

OISTO* 

RESIS. 

AJV (26) • 

PLATE844 

1 

Y* 

R2. 

CT. 

INCY. 

NO 1ST 2. 

INCXO. 

KSENS*. 

AJH(26) * 

PLATE345 

2 

K* 

R3* 

alt  • 

AINCX. 

NDIST3. 

TSTOP. 

nsetup* 

V (26) * 

PLATE346 

3 

T* 

R4» 

NT. 

AINCY. 

N0IST4. 

TEDIT. 

NSETuP2» 

C(26) • 

PLATE347 

4 

AK* 

RAO. 

ALN. 

INCX 1 . 

OISTl. 

OELED. 

NSETUP3, 

VI (26) * 

PLATE848 

S 

OELT. 

ALl. 

AL3* 

INCYI » 

0IST2* 

TE01T2. 

NSETUP4, 

AL (26) . 

PLATE849 

6 

KSEN. 

AL2. 

ANU. 

UELX  • 

0IST3. 

NPOSXl, 

NSETUP5* 

•A (26,26) . 

PLATE850 

7 

NCAP* 

NIOO. 

N2AP. 

OELY. 

0IST4* 

NP0SX2. 

IIIUO). 

•6(26.26)  . 

PLATE351 

8 

XI ( 100) 

• YK100), 

AT(l).  < 

CTRAC(I) 

• CTRK1) 

. CTR2U) 

PLATE852 

01  MENS  ION  AJV  J ( 26*26)  . AJH0(26.26>*  VQ(26.26).  C0(26*26). 

1 AL0QI26.26) 

LEVEL  3.  AJVQ,  AjnQ.  VO,  CQt  ALQO 

COMMON  /SCSI/  AJVQ 
COMMON  /ECS2/  AJrlQ 
COMMON  /ECS 3/  VO 
COMMON  /ECS-*/  CJ 
COMMON  /ECS 5/  AL JO 

COMMON  HAOIUSU) 

PI*3. 14159265 
XMN*YMN*0.0 
XMX*YMX*AMAX1 (X.Y) 

IP(KLM.EO.l)  GO  TO  3 

CALL  LINPLOT(XI.YI.A.XMX.XMN.YMX.YMN) 

return 

CONTINUE 

CALL  GBIDQ  (XMX.XMN.YMX.YMN.XMX.XMN.YMX.YMN* 1 0* 10* 1 • 0. 1HX . 1 . 1HY * 

RADIUS  1 1) aRl 

RADIUS  12) «R2 

radius  <3) *R3 

RAOIUS(A) »rao 

N*50 

DO  2 1*1*4 

ThETA*0.0  

ir  (RADIUS! I) *£Q*0*>  GO  TO  2 
DO  i J*l • V 

XKJ)*x-RaOIuS(I>*SIN(THETa> 

Y| ( J) *Y— RADIUS ( I) »COS (TMETA) 

TMETA*TH£TA«PI/ C2.«N> 

CONTINUE 

CALL  LINPLOT  (XI*Yl.N.XMX.XMN.VMX*YMN) 

CONTINUE 

RETURN 

ENO 


PLATE853 
PLATE854 
PLATE85S 
PLATE3S6 
PLATE357 
PLATE858 
PLATE859 
PLATE360 
PLATE861 
PLATE862 
PLATE863 
PLATE864 
PLATE 865 
PLATE866 
PLATE867 
PLATE36d 
PLATE869 
PLATE870 
PLATE371 
PLATE872 
DPLATE873 
PLATE374 
PLATE875 
PLATE876 
PLATE377 
PLATE373 
PLATE879 
PLATE880 
PLATE881 
PLATE382 
PLATE883 
PLATE884 
PLATE385 
PLATE886 
PLAIE387 
PLATE888 
PLATE889 
PLATE390 
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✓ 


SUBROUTINE  PRINT1  (NEDIT) 


COMMON  X* 

Rl* 

VO* 

INCX. 

NOISTl* 

OISTO* 

RESIS. 

Ajv (26) • 

I 

Y* 

R2» 

CT. 

INCY. 

N0IST2, 

INCXO* 

KSENSa* 

AJH(26> • 

2 

K* 

R3* 

ACT  • 

AINCX* 

NO 1ST 3* 

TSTOP. 

nsetup. 

V (26)  * 

3 

T* 

RA* 

NT. 

AINCY. 

NOISTa* 

TEOIT* 

NSETUP2. 

C (26) • 

A 

AK« 

RAO* 

ALN* 

INCX 1 • 

OISTI. 

OELEO* 

NSETUP3. 

VI  (26)  * 

5 

OELT, 

AL1  • 

ALO* 

INCYI* 

0IST2. 

TE0IT2. 

NSETUPA. 

AL (26) . 

6 

KSEN* 

AL2. 

ANU* 

UELX* 

0IST3* 

NP0SX1* 

NSETUPS* 

•A (26.26) • 

NCAP. 

N100* 

NZAP. 

OELY. 

oista* 

NP0SX2* 

111(10). 

■8(26.26) * 

XI  ( 100) 

• YMIOO). 

AT  ( 1 > . 

CTRACd) 

* CTRK1 

)*  CTR2 ( 1 ) 

01 HENS  I ON  AJV0(26.26>.  AJhO(26.26>*  V0I26.26)*  CJ (26*26). 
1 ALOO (26*26) 


LEVEL  3*  AJVO*  AJHQ*  VO.  CO*  ALOO 

COMMON  /ECS  1/  AJVO 
COMMON  /ECS 2/  AJHO 
COMMON  /ECS 3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECS5/  ALOO 


t 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


t 

2 


This  SUBROUTINE  PRINTS  VARIOUS  EOITEO  ARRAYS  DEPENDING  ON  the 
CALLEO  variable  NEOIT. 

This  SUBROUTINE  MUST  BE  CALLED  »ITh  NEOIT  SET  TO  AN  INTEGER  1 
Through  5. 

AMEN  NEOIT  >1  THE  VOLTAGE  ARRAY  IS  PRINTEO. 

•HEN  NEDIT  «2  The  CAPACITANCE  ARRAY  IS  PRINTEO. 

•MEN  NEOIT  *3  THE  INOUCTANCE  ARRAY  IS  PRINTEO. 

•HEN  NEOIT  •*  Th£  VERTICAL  CURRENT  ARRAY  IS  PRINTED. 

•HEN  NEOIT  as  THE  HORIZONTAL  CURRENT  ARRAY  IS  PRINTEO. 


00 

00 


PRINT  6 
NYa( INCY* 
I r INEOIT 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
2 I ■ I . 

1 Jalt 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
IF  (NEOIT 
CONTINUE 
PRINT  7, 
CONTINUE 
taINCXI 

00  3 jal. 


1*) /IS 

.ES.l)  PRINT  13*  INCX*INCY.INCX1«INCYI 

•CO. I)  PRINT  8 

.E0.2)  PRINT  9 

.EG. 3)  PRINT  10 

.EG. A)  PRINT  II 

•EG.S)  PRINT  12 

INCXltS 

INCY1 

.EG. I)  CALL  REAOEC  (V l ( J) . VO ( I . J) • l) 
.EG. 2)  CALL  REAOEC  (VI ( J) .CG ( I . J) . 1 ) 
.EG. 3)  CALL  REAOEC  (VI ( J) . ALOO ( I . J) . I) 
.EG. •>  CALL  REAOEC  ( VI I J> • AJVO ( I . J) . 1 ) 
.EG.S)  CALL  REAOEC  (VI (J) »AJhg ( I* j) * 1) 

( VI ( J) * Jal , INCY ,NY) * VI ( INCY 1 ) 


INCY  l 
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PLATE891 
PLATE892 
PLATES93 
PLATE89* 
PLATE895 
PLATE896 
PLATE897 
PLATE398 
PLATE899 
PLATE900 
PLATE901 
PLATE902 
PLATE903 
PLATE90* 
PLATE905 
PLATE906 
PLATE907 
PLATE908 
PLATE909 
PLATE910 
PLATE91I 
PLATE912 
PLATE913 
PLATE91* 
PLATE9I5 
PLATE916 
PLATE917 
PLATE918 
PLATE919 
PLATE920 
PLATE921 
PLATE92 2 
PLATE923 
PLATE92A 
PLATE925 
PLATE926 
PLATE927 
PLATE928 
PLATE929 
PLATE930 
PLATE931 
PLATE932 
PLATE933 
PLATE93A 
PLATE935 
PLATE936 
PLATE937 
PLATE933 
PLATE939 
PLATE9A0 
PLATE9aI 
PLATE9A2 
PLATE9*3 
PLATE9AA 


if  <Ncoir.ca.il  call  readcc  <viij>»vq<i»j>»h 

PLATE945 

IF  (NC0IT.EQ.2)  CALL  REAOCC  <V1 <J> .CQ< ItJI <11 

PLATE946 

IF  (NCOIT.EO.3)  CALL  SCAOCC  <VKJI<ALOO<1<JI<D 

PLATE947 

IF  <NC0IT.C0.4)  CALL  SCAOCC  < VI  U>  .AJVQ  < I . J)  < 1 > 

PLATE94d 

IF  <NC0IT.C3.5l  CALL  SCAOCC  < V l < Jl < AJrtQ ( I < J) < D 

PLATE949 

3 

CONTINUE 

PLATE950 

PRINT  7.  IVlIJItJ-l.INCT.NYI.WKlNCYl) 

PLATE9SI 

PRINT  5 

PLATE952 

NMINX>MINO <16<INCY1<INCXD 

PLATE953 

00  4 J»l<MNlNX 

PLATE954 

IF  <NCOIT.CQ.U  CALL  SCAOCC  <V  l < 1 1 < VO  < 1 < J)  <NMlNX) 

PLATC955 

IF  INC0IT.E0.2)  CALL  REAOCC  <VKD<C0<1<JI<NHINXI 

PLATE956 

IF  <NC0IT.CQ.3I  CALL  SCAOCC  <VKD<AL30<1<JI  <NHINX) 

PLATE9S7 

IF  <NCDIT.CQ.4)  CALL  REAOCC  <VKD<AJV0<1<JI<NMINXI 

PLATE95B 

IF  <NCD1T.EQ.5)  CALL  SCAOCC  <VKD<AJH0<1<JI<NN1NXI 

PLATE9S9 

PRINT  7<  <V1 < D < I«1<NHINX) 

PLATC960 

* 

CONTINUE 

PLATE961 

RETURN 

PLATE962 

c 

PLATE963 

C 

PLATE964 

s 

FORMAT  <1X</<10X< 15HL0A0  AREA  ARRAY</) 

PLATE965 

6 

FORMAT  (lHl) 

PLATC966 

7 

FORMAT  <1X<16C3.D 

PLATE967 

a 

format  < i ox < i jhvolt age  array</> 

PLATC968 

9 

format  ( 10X< 17hCAPAC1TANCC  ARRAY</) 

PLATC969 

10 

FORMAT  < 10X< 16HIN0UCTANCE  ARRAY</I 

PLATE970 

li 

format  < iox<22hvcrt i cal  current  akray./i 

PLATC9TI 

12 

FORMAT  I10X<24hh0RIZ0NTAL  CURRENT  ARRAY </) 

PLATC972 

13 

FORMAT  <2X<4IS> 

PLATE973 

end 

PLATE974 

OUU UUU 


SUBROUT INC  PRINT2 


PLATE975 

PLATE976 


COMMON  X* 

Rl* 

VO* 

INCX. 

NO I ST  1 * 

OISTO* 

RESTS. 

AJV (26)* 

PLATE977 

1 

Y. 

R2. 

CT. 

INCY. 

NOIST2, 

INCXO* 

xSENSM* 

AJM ( 26) . 

PLATE978 

2 

K* 

R3. 

ALT* 

AINCX. 

NOIST3* 

TSTOP. 

NSETUP. 

V (26)  • 

PLATE979 

3 

T. 

R*. 

NT. 

AINCY. 

NOIST*. 

TEDIT. 

NSETUP2. 

C (26) • 

PLATE980 

6 

AK* 

RAO* 

ALN* 

INCX  1 * 

OIST1. 

DELED* 

NSETUP3. 

VI (26) . 

PLATE981 

5 

OELT. 

AL1* 

ALQ. 

INCY 1 * 

OIST2* 

TEOIT2. 

NSETUP*. 

AL (26) « 

PLATE982 

6 

KSEN* 

AL2* 

ANU* 

UELX. 

DIST3* 

NPOSX 1 • 

NSETUPS. 

MAI  26. 26) . 

PLATE983 

7 

NCAP* 

N100* 

NZAP* 

OELY. 

OIST*. 

NPOSX2* 

111(10). 

MB (26*26) • 

PLATE98* 

8 

XI (100) 

• Y1 (100) * 

AT  ( 1 ) » < 

CTRAC(l) 

« CTRHl) 

♦ CTR2I1) 

PLATE98S 

01  MENS  I ON  AJVO (26*26)  » AJhO(26*26>*  V<3(26«26>.  CJ(26*26>. 
1 ALQQ (26*26) 

LEVEL  3*  AJVO*  A jHQf  VO*  CO*  ALSO 

COMMON  /CCS1/  AJVO 
COMMON  /ECS2/  AJhO 
COMMON  /ECS 3/  VQ 
COMMON  /ECS*/  CO 
COMMON  /ECSS/  ALOO 


This  SUBROUTINE  PRINTS  SEVCRLY  EDITED  voltage*  HORIZONTAL  current 
ANO  VERTICAL  CURRENT  ARRAYS* 

FORMAT  STATEMENTS 

NZ«(INCY»3)/9 

INNZ«INCT1-NZ“2 

NZ>MAX0(NZ*2) 

PRINT  7 
PRINT  8.-T 
PRINT  9 

NN*«  _ 

00  1 I ■ I * INNZ  *NZ 

NNaNN* I 

XI (NN> «I 

CONTINUE 

NNaNN* 1 

Xl(NN)«INCXl 

PRINT  S*  (X I ( I ) • I ■ I *NN) 

PRINT  6 

00  2 JM.10 

JJ»(J-l)*NZ*l 

I r (J.E0.10)  JJaINCYl 

CALL  REAOCC  (AJV  ( 1 1 ».AJVO  ( 1 1 JJ)  . INCX 1 ) 

PRINT  S*  (AJV  (Dtl*!*  INNZ.NZ)  • AJV  ( INCX  I ) 

CONTINUE 
PRINT  6 
PRINT  IS 
00  3 J-l.10 
J J* ( J“ I ) #nZ ♦ l 


PUATE986 
PLATE987 
PLATE988 
PLATE989 
PLATE990 
PLATE991 
PLATE992 
PLATE993 
PLATE99* 
PLATC99S 
PLATE996 
PLATE997 
PLATE996 
.PLATE999 
PLAT  1000 
PLAT  1 001 
PLAT  1002 
PLAT  1003 
PLAT  100* 

plat  1005 

PLAT  1006 
PLAT1007 
PLAT  1006 
PLAT  1009 
PLAT1010 
PLAT  10 1 1 
PLAT10I2 
PLAT1013 
PLAT101* 
PLAT101S 
PLAT  1016 
PLAT  1017 
PLAT  1018 
PLAT1019 
PLAT  1020 
PLAT1021 
PLAT  1022 
PLAT  1023 
PLAT102* 
PLAT  1025 
PLAT  1026 
PLAT  1027 
PLAT1026 
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p 

* 


i r u.EQ.io)  jj«incyi 

CALL  READEC  ( AJm(1 ) , AJMO< 1 , JJ> • INCx 1 I 

PRINT  St  (AJH(I)  .1«1»INNZ.NZ>  .AJHUNCX1) 

3 CONT INU6 

PRINT  6 
PRINT  11 
DO  A J»1 » 10 
JJ«<J-1)#NZM 
ir  u.EQ.io)  jj«iNcn 
CALL  REAOEC  <V<1).VQ«1.JJ).INCX1) 

PRINT  5.  I V 1 1 ) . !■  1 . INNZ.NZ) t V ( INCA  1 ) 

A CONTINUE 

PRINT  6 
RETURN 

5 FORMAT  I 1*. IP1 1E12.2) 

6 FORMAT  <1N0I 

7 FORMAT  < IhI) 

0 FORMAT  ||0A,24mTI«E  AFTER  INITIATION  ■ .1PE20.8./) 

9 FORMAT  (10A.30KPARTIAL  VERTICAL  CURRENT  ARRAY */) 

10  FORMAT  1 10A. 32NPARTI AL  HORIZONTAL  CURRENT  ARRAY./) 

11  FORMAT  C l OX. 2 IhP ART  I AL  VOLTAGE  ARRAY./) 

ENO 


PLAT1029 

PLAT1030 

plat  1031 
plat  1032 

PLAT1033 
PLAT103A 
PLAT  1035 
plat 103p 
PLAT  1037 
PLAT  1036 
plat  1039 

PLAT1043 
PLAT  1041 
PLAT  1042 
PLAT  1043 
PLAT  104*. 
PLAT1045 
PLAT  1046 
PLAT  1047 
PLAT  1048 
PLAT1049 
PLAT  1050 
PLAT  1051 
PLAT  1052 


V 


SUBROUTINE  INDUCT 


PLAT  1053 
PLAT105A 


COMMON  X* 

Rlt 

VO* 

INCX* 

NO  I ST  1 . 

OISTO* 

RCSIS, 

AJV (26) « 

PLAT  1055 

Y « 

R2. 

CT. 

INCY. 

N0IST2, 

1NCX0. 

KSENSW, 

AJM(26) * 

PLAT  1056 

K* 

R3t 

ALT* 

AINCX* 

N01ST3* 

TSTOP* 

nsctup. 

V (261 « 

plat  1057 

T. 

R4, 

NT* 

A INCY. 

NOISTa* 

TEOIT. 

NSETUP2, 

C (26) • 

PLAT  1058 

AKt 

RAD. 

ALN* 

INCX]. 

DISTl, 

DELED. 

NSETUP3. 

VI (26) « 

PLAT  1059 

OCLT, 

AL1. 

ALQ* 

INCY 1 * 

0IST2. 

TC0IT2, 

NStTUPA, 

AL (26) • 

PLAT  1060 

KSEN. 

AL2. 

ANU* 

UCLX* 

0IST3* 

NP0SX1. 

NSETUP5, 

•A (26*26) • 

PLAT  1061 

NCAP. 

‘(100  • 

NZAP* 

UELY. 

0IST4. 

NP0SX2* 

I II ( 10) * 

WB (26*26) • 

PLAT  1 062 

*1(109) 

• Y 1 ( 1 00 ) • AT ( 1)  ♦ 

CTRAC(l) 

• CTRK1) 

. CTR2 ( 1) 

PLAT1063 

PLAT106A 

DIMENSION 

AJVC (26.26) 

• AJhQ (26*26). 

VQ (26*26) . CQ (26* 

26). 

PLAT  1065 

1 ALOO (26*26) 

LEVEL  3.  A JVC*  AJMQ*  VQ.  CO*  ALOO 

COMMON  /SCSI/  AJVO 
COMMON  /CCS2/  AJMQ 
COMMON  /ECS 3/  VO 
COMMON  /ECSa/  CO 
COMMON  /ECSS/  ALQO 


This  SUBROUTINE  sums  The  VERTICAL  CURRENT  IN  TmE  plate  at  a 
POSITION  A-OISTO. 

IPIT.LE.9.SC-09)  KSEN-0 
KSENaKSEN*! 

KSCNSWaKSENSV*! 

If  (KSCNSW.EQ.1I  PRINT  8 

IP  (KSENSW.CQ.il  PRINT  9 

laOISTO/OELA 

CUR1-O.0 

00  i Jal.INCYl 

CALL  REAOCC  ( A JV ( I ) * AJVQ (I.JItl) 

If  (J.EQ.l)  AJV(I)»AJV(I 1/2.8 

CUR I «CUR 1 • AJV ( I I 

CONTINUE 

If  (CuRl.LC.l.OC-99)  00  TO  2 

ALL»V0*T/CUR1 

CONTINUE 

This  SUBROUTINE  SUMS  ANO  PRINTS  The  CURRENT  CUMING  ' UT  Of  THE 
CAPACITOR  BANKS. 

CURaCUR A>CURB«CURC*9. 9 
If  (NSCTUP.NC.SI  00  TO  « 
laINCXI 
00  3 Ja|.INCY 

CALL  REAOCC  (AJV (I I » AJVQ( I » J> • 1 I 
If  (J.EQ.ll  AJV ( I ) a AJV ( 11/2.9 
CURAaCURA’AJV(l) 

CONTINUE 


PLAT  1066 
PLAT  1067 
PLAT  1066 
plat  1069 
PLAT1070 
PLAT1071 
PLAT1072 
PLAT  1073 
plat  1074 
PLAT107S 
PLAT  1076 
PLAT  1 077 
PLAT 1076 
BLAT1079 
PLAT  1 080 
PLAT1081 
PLAT1082 
PLAT  1083 
PLAT  1084 
PLAT  1085 
PLAT  1086 
PLAT  1 087 
PLAT  1086 
PLAT1089 
PLAT  1090 
PLAT1091 
PLAT1092 
PLAT1093 
PLAT1094 
PLAT  1095 
PL AT  1 096 
PLAT  1 09 ) 
PLAT  1098 

plat  1099 
plat  1100 

PLA t 1101 
PLAT  11 02 
PLAT1103 
PLAT 1104 
PLAT110S 
plat  1 1 0 6 


m uuuuu  r»  uuuu  uu«» 


NN-N01ST1 

PLAT 

NNN«N0IST2 

PLAT 

J-INCY1 

PLAT 

CALL  RCA0CC  ( AJH ( 1 ) . AJHQ ( 1 * J) . INCX) 

PLAT 

DO  A laNN.NNN 

PLAT 

CURBaCURB*AJH (1) 

PLAT 

CONTINUE 

PLAT 

NNaNOISTO 

PLAT 

NNNaNOIST  A 

PLAT 

OO  5 laNN.NNN 

PLAT 

CURCaCURC*AJH(i> 

PLAT 

CONTINUE 

PLAT 

CONTINUE 

PLAT 

CU9aCURA.CUR8.CJRC 

PLAT 

PLAT 

this  subroutine  sums  the  charge  across  each 

CAPACITOR  ANO  PRINTS 

plat 

out  the  result  as  total  charge,  this  allots 

CONSERVATION  OF 

PLAT 

charge  to  6E  checked. 

PLAT 

PLAT 

CHaO.O 

PLAT 

00  7 Jal.INCYl 

PLAT 

CALL  REAOEC  (Vll).VO(l.J). INCX 1 ) 

PLAT 

CALL  REAOEC  (C « 1 1 .CO ( 1 . J) . INCX1 1 

PLAT 

00  7 lal.INCXl 

PLAT 

IF  (J.GT.I)  GO  TO  7 

PLAT 

IF  (J.EO.l)  V<I>aV(I»/2.0 

PLAT 

CH*CH.V ( I ) *C ( I ) 

PLAT 

CONTINUE 

PLAT 

PRINT  10.  T.  CUR 1 « CUR .CH. ALL ‘CURA. CURB. CURC 

PLAT 

PLAT 

This  SECTION  SAVES  CURRENTS  AND  TIMES  FOR  A 

time  oepenoent 

PLAT 

CURRENT  TRACE  TO  BE  PLOTTCO  BY  SUBROUTINE  CTRACE. 

PLAT 

.. 

PLAT 

AT (KSEN)aT 

PLAT 

CTRAC (KSEn) »CUR 1 

PLAT 

CALL  INDUCT 2 (CURD 

PLAT 

return 

PLAT 



PLAT 

PLAT 

format  (IhD 

PLAT 

FORMAT  (10X.AHTIME.BX.AHCURl.8X.AHCUR2.7X.6nCHARGE.6X.6HlNOUCT.7X.PLAT 

1AHCURA • 8* . OHCUPB t ax . bnCURC . / ) 
10  FORMAT  (5X.IP8C12.ai 
CNO 


PLAT 

PLAT 

PLAT 


1107 

1108 

1109 

1110 
1111 
1112 
1113 

1 1 1 A 
1115 
111b 

1117 

1118 

1119 

1120 
1121 
1122 
1123 
112a 
1125 
112b 

1127 

1128 

1129 

1130 

1131 

1132 

1133 
113a 

1135 

1136 

1137 

1138 

1139 
1 1A0 
1141 
11A2 
11  A3 
IU* 
HAS 
11  AO 
11A7 
11A8 
11A9 

ISO 
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i 


SUBROUTINE  INOUCT2  (CURD 


PLAT1151 
PLAT  1152 


COMMON  A* 

• 1* 

VO* 

INCA* 

NO  1ST  1 * 

OISTO* 

RESIS. 

AJVI26) « 

PLAT1153 

1 

Y* 

R2. 

CT. 

INCY. 

N0IST2* 

INCA9* 

KSENSa* 

AJH(2b) « 

PLATUS* 

2 

K* 

R3. 

ALT. 

A INC  A* 

N0IST3* 

TSTOR. 

NSETUP, 

V (26) « 

PLAT  1 155 

T. 

R6* 

NT* 

AINCY. 

N0IST4. 

TE01T* 

NSETUP2* 

C (26)  • 

PLAT11S6 

A*. 

RAO* 

ALN* 

INCA  1 • 

0IST1. 

OELEO* 

NSETUP3. 

VI (26) • 

PLAT  1157 

XLT. 

AL1* 

ALO* 

INCYl. 

0IST2. 

TE0IT2. 

NSETUP*. 

AH 26)  • 

PLAT  1158 

• SEN* 

AL2* 

ANU* 

OELA* 

0IST3, 

NPOSAl, 

NSETUP5. 

•A(26«26> • 

PLAT  1159 

NCAP* 

N100* 

NZAP* 

UELY  • 

OIST*. 

NP0SA2* 

111(10)* 

•6(26*26)  * 

PLATU60 

XI (104) 

« Y I ( 1 'JO  > • 

AT  ( 1 > • 

CTRAC ( I) 

. CTRKl) 

. CTR2 ( 1) 

PLATU61 

01 MENS I ON  AJVO (26.26) • AJhU(26*2o)*  73(26*26)*  03(26*261* 
1 ALOC  (26.26) 

LEVEL  3*  AJVO*  AJH3.  VO.  CO*  ALOQ 

COMMON  /ECSl/  AJVO 
COMMON  /ECS2/  AJHO 
COMMON  /ECS3/  VU 
COMMON  /ECS*/  CO 
COMMON  /ECS5/  ALSO 

M10«M10*1 

M10«MOO(M10.10) 

if  (Mio.NE.oi  return 

VOLTI-VOLT2 

TIM1.TIM2 

CU1*CU2 

CALL  RCAOEC  ( V0LT2* VO ( INCX1 **) « 1) 

T JM2«T 
CU2-CUR1 

IF  (CU2.E3.CUP  RETURN 

ALNEv» ( ( VOL  T 1 * V0LT21 /2) • ( T I M2-T I M 1 1 / ( CU2-CU 1 1 
PRINT  I*  ALNEV 

RETURN  - 


PLAT  11 62 
PLAT1163 
PLAT  116* 
PLAT116S 
PLATU66 
PLAT1167 
PLAT  1166 
PLAT1169 
PLAT 11 73 
PLAT1171 
PLAT  11 72 
PLAT  11 73 
PLAT117* 
PLAT117S 
PLAT  11 76 
PLAT  1177 
PLAT  1178 
PLAT  11 79 
PLATliaO 
plat liai 
PLAT1182 
PLAT1183 
PLAT  118* 
PLAT l 135 
plat l 186 


PLATU87 

FORMAT  ( 10Xr62HTME  AVERAGE  INOUCTANCE  SINCE  THE  LAST  STATEMENT  LIKPLAT 1183 
IE  This  is  .E12.6.10H  HENRIES.*/)  PLATU89 

eno  plat  1190 


SUBROUTINE  CTRACE 


PLAT1 191 
PLAT  1192 


COMMON  X* 

Rl* 

vo« 

INCX* 

NO  I ST  1 * 

OISTO* 

RESIS. 

AJV (261  * 

PLAT  1193 

1 

Y* 

R2* 

CT* 

INCY. 

M0IST2* 

INCXO* 

KSENSv* 

AJH(26>  * 

PLAT  1194 

2 

K. 

R3» 

ALT. 

AINCX* 

NQIST3* 

TSTOP. 

nsetup. 

V (261  * 

PLAT  11 95 

3 

T* 

R4, 

NT. 

AINCY* 

N0IST4* 

TE01T* 

NSETUP2. 

C (26) < 

PLAT  1196 

A 

AK* 

RAO* 

ALN* 

INCX 1 • 

0IST1, 

OELEO* 

NSETUP3* 

VI  (26)  * 

PLAT  1197 

s 

OELT* 

AL1* 

ALO* 

INCY 1 * 

0IST2. 

TE0IT2* 

NSETUP A, 

AL  ( 261  « 

PLAT  1198 

6 

KSEN* 

AL2* 

ANU* 

OELX* 

0IST3. 

NP0SX1* 

NSETJP5* 

■A (26*26) * 

PLAT  1199 

7 

NCAP* 

N100* 

NZAP* 

OELY* 

OISTA* 

NP0SX2* 

iimoi* 

■8(26*26) * 

PLAT  1200 

8 

XI (1001 

• YK100I* 

AT  ( 1)  * 

CTRACtll 

• CTRK1 

1*  CTR2(11 

Plat i2o i 

DIMENSION  AJVQ (26*261*  AJhQ (26*261*  VO (26*261*  CO (26* 
1 ALOO (26.261 

LEVEL  3*  AJVQ*  AJHQ*  VO*  CO*  ALOO 

COMMON  /ECS1/  AJVQ 
COMMON  /ECS2/  AJHQ 
COMMON  /ECS 3/  VO 
COMMON  /ECS**/  Cu 
COMMON  /ECSS/  ALOO 

XMN»XMX»YMN»YMX«0.0 
IF(KSEN.LT.3)  RETURN 
00  1 K«1»KSEN 
XMN>AMINI(XMN<AT(K11 
VMNaAMINl (YMN.CTRAC (K1 1 
XMXaAMAX 1 ( XMX . AT (K) ) 

YMX»AMAX1 (YMX.CTRAC (K1 I 
CONTINUE 

CALL  GRIOQ  (XMX*XMN*YMX*YMN*XMX*XMN*YMX*YMN*10*10«1*0* 
1URRENT  * 71 

CALL  LINPLOT  < AT.CTRAC.KSCN.XmX. XMN.YMX. YMN) 

CALL  LINPLOT ( AT *CTRaC»KSEN* amx.XMN* ymx.YMN) 

CALL  FRAME 
RETURN 
END 


PLAT  1202 
261*  PLAT  1203 

PLAT  120** 
PLAT120S 
PLAT  1206 
PLAT1207 
PLAT1206 
PLAT1209 

plat  12 1 o 

PLAT1211 
PLAT1212 
PLAT1213 
PLAT1214 
PLAT1215 
PLAT1216 
PLAT  121 7 

plat  1218 

PLAT1219 
PLAT1223 
PLAT  1221 
AMTIME.A. 7MCPLAT 1222 
PLAT  1223 
PLAT122** 
PLAT122S 
PLAT  1226 
PLAT  1227 

plat  1228 
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SUBROUTINE  CTRAC2 


1 


COMMON  X* 

r. 

Kf 

T. 
AK* 
OELT  * 
KSEN* 


Rlt 

R2* 

R3* 

R4. 

Rad* 
AL1  * 
AL2* 


VO* 

CT* 

ALT* 

NT* 

ALN* 

ALO* 

ANU* 


NCAP*  NlOO*  N2AP* 


INCX*  NOISTl*  OISTO*  RESIS, 

INCY*  N0IST2.  INCXO*  KSENSV* 

AINCX*  N0IST3*  TSTOP.  NSETUP* 

AINCY.  NOISTa,  TEOIT.  NSETUP 2. 

INCX  I • OISTl*  OELEO*  NSET0P3* 

INCY l*  0IST2,  TE0IT2.  NSETUP4, 

OE LX*  0IST3*  nPOSX I * NSETUP5*  «A(26.26). 
OELY*  OISTA*  NPOSX 2*  1 1 1 ( 10> • 46(26*26)* 


PLAT  1229 
PLAT  1230 
AJV (26) * PLAT1231 
A JH ( 26) « PLAT  1232 
V (26) « PLAT  1233 
C (26) * PLAT  123a 
V 1 1 26) • PLAT123S 
AL (26) • PLAT1236 
PLAT  1237 


Xl(lOO)*  Y1  ( 190)  • AT  ( 1)  • CTRACO)  * CTRl(l)*  CTR2I1) 


DIMENSION  AJVQ (26*26) • AJhO(26.26).  V0(26*26)*  CJ(26*26>. 
1 ALQO (26.26) 

LEVEL  3*  AJVQ*  AJmO*  VU.  CO*  ALOO 

COMMON  /ECS1/  AJVQ 
COMMON  /ECS2/  AJMQ 
COMMON  /ECS3/  VQ 
COMMON  /ECSa/  CO 
COMMON  /ECSS/  ALOO 


XMX«XMN»YMN«YMX»0.9 
IFIN190.lt. 3)  RETURN 
DO  l *«1.N100 
XMNaiM INI (XMN.AT(K) ) 

YMN»AM1N1 (YMH.CTRl (K) ) 

XMXAAMAX1 (XMX • AT (K ) > 

TMX«AMAX1 (TMX.CTR1 (K>) 

CONTINUE 

TMA»1.S»YMX 

Call  or i oo  (xmx*xmn»vma*ymn*xmx.xmn.ymx.ymn.io* io*i*o*ahtime.a,7hcplat1261 

1URM001.7)  - PLAT  1262 

CALL  LINPLOT  <AT*CTR1*N10«.XMX.XMN,YMX,YMN)  PLAT1263 

CALL  LINPLOT ( AT  «CTR l *N100  ,XMX*XMN,YMX.YMN)  PLAT126A 

CALL  FRAME  PLAT1265 

CALL  OR  I DO  (-AMX.XMN.  YMX  * YMN*  XMX  * XMN.  YMX*  YMN*  19*10*1*  0*  AMT  IM£*  A,  7HCPLAT 1266 


PLAT  1238 
PLAT  1239 
PLAT12A0 
PL AT  12a 1 
PLAT12A2 
PLAT12A3 
PLA71244 
PLAT12AS 
PLAT12A6 
PLAT 12A 7 
PLAT1248 
PLAT  1249 
PLAT12S0 
PLAT  1251 
PLAT12S2 
PLAT1253 
PLAT  1 25- 
PLAT  1255 
PLAT  1256 
PLAT1257 
PLAT12S8 
PLAT12S9 
PLAT  1260 


1URM002*  7) 

CALL  LINPLOT  (AT*CTR2*N100*XMX*XMN* VMX.YMN) 
CALL  LINPLOT (AT  *CTR2*N109  « XMX. XMN. YMX. YMN) 
CALL  FRAME 

return 

ENO 


plat  1257 
PLAT  1266 
PLAT  1269 
PLAT1270 
PLAT1271 
PLAT  1272 
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c 

c 

c 

c 


subroutine  current 

DIMENSION  LABEL (2) • XTRALB (2) 


PLAT  1273 
PLAT1274 
PLAT1275 


COMMON  X* 

Rlt 

VO# 

INCX# 

N0IST1# 

01STO# 

RESIS, 

AJVI26) • 

plat  1276 

1 

Y. 

R2» 

CT. 

INCY# 

NO  I S T 2 * 

INCXO# 

KSENSa. 

AJH(26> • 

plat  1277 

2 

K. 

R3* 

ALT# 

A INCX# 

N0IST3. 

TSTOP. 

NSETUP. 

V (26) • 

plat  1 278 

3 

T. 

AM. 

NT# 

AINCY. 

NOISTM. 

TEDIT. 

NSETUP2. 

C (26) • 

PLAT1279 

4 

AK# 

RAO# 

ALN# 

INCX 1 * 

DIST 1 • 

OELEO# 

NSETUP3# 

V 1 (26) • 

PLAT  1230 

5 

OELT. 

AL1# 

ALO# 

INCY 1 • 

0IST2. 

TE0IT2. 

NSETUPm, 

AL (26) # 

PLAT1281 

6 

KSENf 

AL2# 

ANU# 

DELX# 

0IST3# 

NP0SX1# 

NSETUP5. 

WAI26.26) . 

PLAT  1282 

7 

NCAP. 

N100# 

NZAP* 

OELY# 

0IST4. 

NP0SX2. 

111(10) # 

Ad(2o#26> . 

PLAT1283 

a 

Xl(lOO)#  Y 1 1 1001 # 

ATI11# 

CTRACUl 

. CTRK1I 

• CTR2M) 

PLAT  128m 

DIMENSION  AJVQ (2S#26>  * AJhU<26.26).  V3(26»26>.  C0(26#26># 
l ALOQ 126.26) 

LEVEL  3*  AJVOt  AJMQ#  VO#  CO#  ALOQ 

COMMON  /ECSl/  AJVO 
COMMON  /ECS2/  AJrtU 
COMMON  /ECS3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECSS/  ALOO 

COMMON  AX (1000)  tATUOOO) 


PLAT  1285 
PLAT  1286 
PLAT1287 
PLAT1288 
PLAT1239 
PLAT  1290 
PLAT  1291 
PLAT1292 
PLAT  1293 
PLAT129M 
PLAT1295 
PLAT  1296 
PLAT1297 
PLAT1298 

this  subroutine  follows  several  current  paths  in  The  transmission  plat i 299 

E0R  A SPECIFIC  time  Th£M  PLOTS  THESE  PATHS  ON  MICROFILM.  PLAT1300 

PL AT  1 301 

XHN»YMN«0.0  PLAT  1302 

YMXaXMXaAMAX 1 ( X t Y ) PLAT1303 

IF  <NSETUP.E0.S1  YMX«XNX*I.5*Y  PLAT1304 

LABEL <1)«5HTIM£»  PLAT  1305 

LABEL <21 «*H  SEC  - PLAT1306 

ENCOOE  (20»7*xTRALB1  LABEL  1 1) tT. LABEL <2)  PLAT1307 

NKT»2  - PLAT  1308 

PI2*3« 161592653S/2.0  - PLAT1309 

IF  (nsetup.ne.5)  naT»i  plat  1310 

oo  6 mmm«i»nkt  plat  131 1 

IF  (MMM.E0.2T-XHXaYMX>X  . PLAT13I2 

AlNCRa  <XMX*YMX) / 1 000.  PLAT1313 

CALL  GRIOO  <XMX.XMN.YMX.YMN#XMA.XMN.YMX.YMN#10.10.1*1#1HX.1.1HY.I.PLAT131m 
1XTRAL8.20)  PLAT1315 

MNKa21  PLAT1316 

IF (NSETUP.EQ. 1)  MNKaINCY 1 PLAT1317 

DO  S L-I.mNK  PLAT  1318 

FIXaMNK-l.  PL AT  13 19 

AXMlaXMX  PLAT  1320 

AT<l»a<L-l)*T/FIX  PLAT1321 

IF  INSETUP. NE, 5)  60  TO  1 PLAT1322 

LLa<L*l»/2  PLAT  1323 

LQ«L/2  PLAT  1329 

IF  IMMM.EO.l.OR.LO.LT.LL)  60  TO  1 PLAT132S 

AT  Ml  ay  PLAT  1326 
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AX(1)«( (L-11/FIX)«X 

PLAT1327 

CONTINUE 

plat  1323 

IF  (NSETUP.NE.l)  60  TO  2 

PLAT  1 329 

ThET»(L-l) *P I 2/FIX 

PLAT1330 

AX ( 1) «6.85#X«C0S(Th£T1 

plat  1331 

AY(I)«0.35*Y*SIN<TnET) 

PLAT1332 

CONTINUE 

plat  1333 

00  3 Ka 1,999 

PL AT  133— 

Ia(AX(K)/X)«(AINCX-U>  *2. 

PLAT133S 

Ja(AY(X)/Y)*(AINCY-l.»  *2. 

PLAT  1 336 

X3»( 1-2) *0£LX 

PLAT  1337 

X2a(I-l)«0ELX 

PLAT  1338 

Y 3* ( J— 2) »OELY 

PLAT  1339 

Y2* ( J- l ) *0£LY 

PLAT  1340 

CALL  PEAOEC  (AJii(I-l)t  A JHO  ( I-l  * J)  *2) 

plat i 3a l 

CALL  REAOEC  <AJV<2>.AJV0<I.J>*1> 

PLAT13A2 

CALL  PEAOEC  (AJV ( 1 ) « AJVQ ( I » J-l ) * 11 

PLAT1343 

CURMbAJH ( I ) • ( X2-AX (K) I /OELX • AJrt ( I- 1 ) * ( AX  (X ) -X3> /OELX 

PLAT  134- 

CURVaAJV (21 • ( Y2-AY (K) ) /UELY* AJV ( 1) • 1 AY (X) -Y3) /CELY 

PLAT13-5 

CUR«SOST (CURH«CURH*CUSv*CURV) 

PLAT  1 346 

IF  (CUR. LE.  1.0)  GO  TO  4 

PLAT  1347 

AX  (X ♦ I)  «AX  (X)  -AINCH*C'JRV/CUR 

PLAT  1346 

AY(X*llaAY(X)-AINCR*CURM/CUR 

plat i 349 

NaX 

PLAT  1350 

IF  UX<K*i)  .LT.O.O.OR.AX(K»1>  .6T.XMX.0R.AY<K*1)  .LT . 0. O.OR . AY (X . 1 > .PLAT  1351 
1ST.YMX)  60  TO  4 PLAT13S2 

CONTINUE  PLAT  1353 

CONTINUE  PLAT  135- 

CALL  LINPLOT  (AX.AY.N*AMX.XMN.YMX.YMN)  PLAT1355 

CALL  LINPLOT<AX*AY*N*XMX»X>«N»YMA*TMN>  PLAT1356 

CONTINUE  PLAT  1357 

CALL  FRAME  PLAT1353 

CONTINUE  PLAT  1359 

RETURN  - PLAT  1360 

PLAT  1361 
PLAT l 362 

FORMAT  (A3* 1PE1 1 .2* A4>  ~ PLAT1363 

END  PLAT  1 36- 
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subroutine  unwind 


PLAT  1365 
PLAT  1366 


COMMON  X. 

RI. 

VO. 

INCX. 

N0IST1, 

OISTO. 

RESIS. 

AJV (26)  . 

PLATI367 

1 

Y. 

R2. 

CT, 

INCY. 

NO  I ST 2* 

INCXO. 

KSENS*. 

AJM<26) . 

PLAT  1368 

2 

K. 

R3. 

ALT. 

AINCX. 

N01ST3. 

TSTOP. 

NSETUP. 

V (26)  • 

PLAT  1369 

3 

T. 

R*. 

NT. 

A INCY. 

NO  I ST*. 

TEOIT. 

NSETUP2. 

C (26) . 

PLAT1370 

* 

AK. 

RAD. 

ALN. 

INCX 1 • 

0IST1, 

OELEO. 

NSE7UP3. 

VI (261 • 

PLA71371 

S 

OCLT. 

AL1. 

ALO. 

INCY 1 . 

0IST2, 

TEDIT2. 

NSETUP*. 

AL (26) • 

PLAT  1372 

6 

KSEN. 

AL2. 

ANU. 

OELX. 

DIST3. 

NP0SX1. 

NSETUP5. 

a A ( 26. 26 ) • 

PLAT1373 

7 

NCAP. 

N100. 

NZAP. 

JELY. 

OIST*. 

NP0SX2 . 

I I I ( 10) • 

*8(26.26) . 

PLAT137* 

8 

XU  100) 

• YK100). 

AT  ( 1 ) , 

CTRAC ( 1 ) 

. CTRl(l) 

. CTR2U) 

PLAT1375 

DIMENSION  AJVO (26.26) » AJMQ(26.26>.  VQ<26.26>.  C3(26.26>. 

I ALQQ (26. 26) 

LEVEL  3*  AJVO.  AJhQ.  VQ»  CO.  ALQQ 

COMMON  /ECS1/  AJVO 
COMMON  /ECS2/  A JHU 
COMMON  /ECS3/  VO 
COMMON  /ECS*/  CO 
COMMON  /ECS5/  ALQQ 

COMMON  OELPmI (509) » CJR(500>«  THETA(SOl).  CURAV(SOO) 

01 MENS  I ON  OELPMIA<500> 

This  SUBROUTINE  COMPUTES  Th£  RELATIVE  CURRENT  DENSITY  CROSSING 
AN  ARC  OUTSIDE  Trt£  LOAO  AND  IT  PLOTS  CURRENT  VERSUS  THETA 
(IN  RAOIANS).  IT  THEN  PLOTS  PH  VESSJS  Th£TA.  MH£9E  PmI  IS  Th£ 
ANGLE  BETWEEN  ThE  CURRENT  VECTOR  AhO  THE  NEGATIVE  RAOIUS  VECTOR. 

EQUIVALENCE  (CURAV 1 1 ) * DELPHI A ( 1 ) ) 

PI 2*3. 1*15926535/2.0 

Th£TA(1>«0.0 

N«S«0 

AN»N 

XMXaPIZ 
AMN»0 . 

ZMX*.50  __ 

Zmn»-zmx 
00  10  LM-1.20 
RAZ«RI»0.02*LM«Y 
00  3 R*l .N 
ay»raz*sin«thcta<k> ) 
ax»raz»cos ( theta <k > ) 
ir  (nsetup.ne.2)  go  to  i 

AVaTHETA(K) 

AX-RAZ 

CONTINUE 

Ia<AA/X)*(AINCX>I.)*2. 

Ja( AV/Y) • > AlNCY» 1 « ) *2. 

A«a<I-2)*CELX 

»*•(/-»)  *«Lf 


PLAT  1376 
PLAT  1377 
plat i 378 
PLAT  1379 
PLAT  1 380 
PLAT  1381 
PLAT  1382 
PLAT  1383 
PL AT  1 38* 
plat  1385 
PLAT  1386 
PLAT1387 
plat i 388 
PLAT13S9 
plat  1390 
PLAT  1391 
PLAT  1392 
PLAT1393 
PL  AT  1 39*. 
PLAT  1395 
PLAT  1 396 
PLAT1397 
PLAT  1398 
PLAT  1 399 
PLAT1*00 
PLATlAOl 
PLAT  1*02 

plat  1*03 

PLAT  1*0* 

plat  1*05 

PLATl*0b 
plat i*o7 
plat i*08 
plat  1*09 
PLAT1*10 
PLAT1*11 
PLAT1*12 
plat i *i 3 
plat i*i* 
PLAT  1*15 
PLATU16 
PLAT  1*17 
PLATU13 


J 


o*  Ul 


* 


2 


3 


4 


7 


Y2»(J-1)*DELY 

CALL  REAOEC  ( AJH ( 1-1 ) . AJMO ( 1-1 . J> *2> 

CALL  REAOEC  ( AJV (2> • AJVU ( I , Jl . 1) 

CALL  REAOEC  (AJV(l),AJVO(J,J-l).l> 

CURHvAJH ( I ) • ( A2-AX ) /OELX*A JH 1 1-1  > * ( AX-XO) /OELX 
CURV«A  JV  ( 2)  * ( Y 2- AY  > /OELY  ♦ A JV  ( 1 ) * t AY-Y  Cl  > /DELY 
CUR (K) "SORT (CURH*CURrt«CURv*CURY) 

IF  (CUR(K).LE,0.)  ?Hl«0. 

IF  <CUR(K)  .LE.O.)  00  TO  2 
PHI*ACOS(CURV/CURIX)  > 

CONTINUE 

DELPHI  «K)»PHI-THETA(lO 
IF  (K.EQ.l)  YMN*YMX*CUR(K> 

Y MN*  AM l N 1 ( YMN  . CUR ( K ) ) 

YMNmfl.O 

YMX*AMAX1(YMX.CUR(K> ) 

THETACK*1>*ThCTA(K> »P 12/AN 

IF  (NSETUP.E3.21  TriETA<<M)*K«Y/AN 

CONTINUE 

IF  (LM.NE.2)  GO  TO  * 

CALL  OR 1 00  ( XMX. XMN* YMX  * YMN* XMX  * XMN. YMX • YMN. 10. 
I CURRENT » 71 

CALL  LINPLOT  (THETA. CUR. N. XMX. XMN. YMX. YMNI 
CALL  LINPLOT (ThETA.CUR.N.XMX.XMN. YMX.YMN) 

CALL  FRAME 

CONTINUE 

CUR  ( 1)  *CUR  ( 2) 

NN*5I 

00  6 K*1.N 
LL»*-l-NN/2 

IF  ILL.LT.-01-  CURAV <K) *CUR (-LL) 

IF  (LL.E0.31  CURAY(K)«CUR(1> 

IF  (LL.GT.O)  CURaV(K)*CuR(LL> 

00  5 1*2. NN 
kk*I*k-2-NN/2 


PL AT  1 A 19 
PLAT  1420 
PLAT1421 
PL AT  1422 
PLAT  1423 
PLAT  1424 
PLAT  1425 

plat i42o 

PLAT1427 
PLAT1428 
PLAT  1429 
PLAT  1430 
PLAT  1431 
PLAT  1432 
PLAT  1433 
PLAT1434 
PLAT143S 
PLAT  1436 
plat  1437 
PLAT  1438 

10. I .0.5HTHETA.5. 7nPLAT 1439 

plat  1440 

PLAT  1441 
PLAT  1442 
PLAT  1443 
PLAT  1 444 
PLAT  1445 
PLAT  1446 

plat  1447 

PLAT  1448 

plat  1449 

plat  1450 
plat  1451 
PLAT  1452 
plat  1453 


KKK*-KK 

IF  (KK.LT.O)  CURSE*CU»(KKK) 

IF  (KK.Ea.O>_CURSE*CUP(l) 

IF  (KK.GT.500)  CU*SE»CUR(2*N-KK> 

IF  (KK.OE.1.ANO.KK.LE.500)  CUR3£“CuR (KK) 
CURAV (<) aCURAV (X) »CURS£ 
CURAV(K)*CURAV(K)/NN 
CURMX  *CUR.MN*CURAV  ( 1 ) 

CURX*CURN*CUR ( 10) 

00  7 K*2.N 

CURMX  * AM AX | ( CUR A V <K ) * CURMX ) 

CURMN*AMIN1 (CURAV (K> .CURMN) 

CURXsAHAX 1 ( CUR ( K ) . CJRX) 

CUPNsAMInI (CUR (X) » CURN ) 

CONTINUE 

OE V I ATE* ( CURMX -CURMN) *200 . / ( CURMX  »CURmN ) 
0EV2* (CURX-CURNI *200./ (CURX*CURN> 

PRINT  11.  R1.R2.RAZ 
PRINT  12.  0EVIATE.0EV2 


PLAT  1*54 
PLAT  1455 
PLAT1456 
PLAT1457 
PLAT  1453 
PLAT1459 
PLAT  14*0 
PLAT1461 
PLAT  1462 
PLAT  1463 
PLAT  1464 
PLAT  1465 
PLAT  1466 
PLAT  1467 
PLAT  1468 

plat  1469 

PLAT1470 
PLAT1471 
PLAT  1472 


p 


80 


IF  ILM.NE.2)  CO  TO  10 

CALL  GRIOQ  (XMX*ANN#YHX*YMN*AMA*AMN* YMXtYMNt 10« 10* 
1CUPAW.S) 

CALL  LINPLOT  ( THET  A«CURAV  tN» AHA • ANNt  YMA* YMN) 

CALL  LINPLOT ( THET  A* CUR AV  *N* AMX  * XMN»  YMX  » YMN) 

CALL  FRAME 

NN»50 

M»N-N N 

00  9 K«1«M 

DELPHI A (K) aOELPHl (K) 

00  8 t«2.NN 

DELPHI A <K ) aOELPHl A IK) *0ELPH1 IK  • 1“  1 ) 

0CLPH1 AIK) «0ELPH I A ( K ) /NN 

CALL  GPIOO  (XMXtA.MN«ZMA*2MN*XMX*XMN»ZMX«ZMN* 10« lOt 
IDELPHI A, 7) 

CALL  LINPLOT  ( Th£T  A t OELPH I A • M t AH A • XHN • ZHX • ZHN ) 

CALL  LINPLOT I THETA. DELPHI  A tM.XMAtA.MN,ZMAiZMN) 

CALL  FRAME 

CALL  GPIOO  IXMX«XMN»ZHX*ZMN.AMA.AHN.ZHX.ZMN.10,10. 
1 DELPHI *61 

CALL  LINPLOT  ( TmET A*OELPHI ,N,X.MX t X.NN*ZMX . ZHN) 

CALL  LINPLOT (THETAtDELPHI  »M t XMX  * AMN*  ZMX  # ZHN) 

CALL  FPAHE 

CONTINUE 

RETURN 


FORMAT  I13A.5"R1  a .E10.3.7H  R2  ■ »E19.3*8h  RAD 
FORMAT  (10A.37mP£RC£NT  THETA  VARIATION  OF  CURRENT 
12  a ,PlOv2»//> 

END 


PLAT  1473 
1 • 0«5HTH£T  At5«  5HPLAT 1474 
PLAT  1475 
PLAT  1476 
PLAT  1477 
PLAT1478 
PLAT  1479 
PLAT1480 
PLAT  1481 
PLAT  1482 
PLAT  1483 
PLAT  1484 
PLAT  1485 
1»0»5HTHETA*5»7mPLAT1486 
PLAT  1487 
PLAT  1488 
PLAT  1489 
PLAT  1490 
1 » 0»5hThET  At 5»6nPL AT  1491 
PLAT  1492 
PLAT  1493 
PLAT  1494 
PLAT149S 
PLAT  1496 
< PLAT  1497 
PLAT  1498 
PLAT 1 499 

■ .Eio.3>  plat  1500 

« .FI0.2.9H  0EVPLAT1501 

PL  AT.  1502 

PLAT1503 


SUBROUTINE  LINPLOT  ( A .d.N. AMAX » AMIN. 3M AX • BH IN) 

PLAT150A 

dimension  Ad),  am 

plat  isos 

c 

PLAT  1506 

c 

THIS  subroutine  plots  data  on  microfilm. 

PLAT  1507 

c 

PLAT  1503 

CALL  SMAPX  fAMIN.AMAX. 0.079. 0.979) 

PLAT  1509 

CALL  SMAPY  (8MIN.3MAX. 0.079. 0.979) 

PLAT  1510 

CALL  SVTRS  ( A.d.N) 

PLAT1511 

return 

PLAT  1512 

ENO 

PLAT1S13 

AFWL-TR-76-205 


APPENDIX  F 
SAMPLE  OUTPUT 

The  Input  for  the  second  example  In  Section  IV  produced  the  printed  out- 
put that  Is  listed  on  the  following  pages.  When  variable  names  are  used,  they 
correspond  to  the  same  variables  In  the  code.  The  definitions  of  most  of  these 
variables  can  be  found  In  the  list  of  variables  In  Appendix  A. 
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